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SUMMARY 


The  final-stage  rocket  which  projected  Skylab  into  orbit  on  14  May  1973 
itself  entered  a  nearly  circular  orbit  at  a  height  near  400  kin,  inclined  at  50 
to  the  equator.  The  rocket,  designated  197 3—2 7 B ,  remained  in  orbit  until 
11  January  1975  and,  being  25  m  long  and  10  m  in  diameter,  was  the  brightest  of 
the  artificial  satellites  then  visible.  Its  orbit  has  been  determined  at  62 
epochs  from  some  5000  optical  and  radar  observations.  The  average  orbital 
accuracy  in  perigee  height  and  orbital  inclination  is  90  m,  though  some  orbits 
with  Hewitt  camera  observations  have  much  smaller  sd  (down  to  10  m) . 

As  the  orbit  contracted  under  the  influence  of  air  drag,  it  passed  slowly 
through  the  31:2  geopotential  resonance,  when  the  track  over  the  Earth  repeats 
every  31  revolutions  at  intervals  of  2  days.  The  variations  in  inclination  and 
eccentricity  during  the  resonance  phase  have  been  successfully  analysed,  after 
some  difficulties,  and  values  have  been  obtained  for  six  lumped  3lst-order 
harmonic  coefficients  in  the  geopotential;  these  will  provide  a  crucial  test  of 
the  accuracy  of  future  geopotential  models. 

The  variation  in  inclination  before  and  after  resonance  has  been  analysed 
to  determine  the  average  atmospheric  rotation  rate  A  (rev/day) .  Results  are: 

A  =  1.04  -  0.05  at  height  380  km  between  June  1973  and  June  1974;  A  =  1.34  -  0.09 
at  height  305  km  for  October-December  1974;  and  A  =  1.06  ±  0.06  at  height  200  km 
in  January  1975.  A  rapid  atmospheric  rotation  at  heights  of  300-350  km,  deduced 
previously  by  combining  results  from  many  satellites,  is  thus  confirmed  by  a 
single  satellite. 

The  observed  variations  ir.  eccentricity  have  been  compared  with  those  pre¬ 
dicted  by  the  (untested)  theory  for  orbit  contraction  in  an  atmosphere  with 
strong  (threefold)  day-to-night  variation  in  density.  3'he  comparison  confirms 
the  accuracy  of  the  theory  and  of  the  CIRA  1972.  atmospheric  model. 

Departmental  Reference:  Space  564 
Copyright  ©  Controller  HMOO London  19? 9 
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INTRODUCTION 
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The  final-stage  rocket  which  fired  Sky  lab  1  inLo  orbit  on  14  May  1471  it; 
one  of  the  largest  objects  ever  to  have  achieved  a  lengthy  orbital  tile.  Desig¬ 
nated  1473-4711,  it  was  a  35-ton  cylinder*  with  u  length  of  35  m  and  a  diameter  ot 
10  m,  so  that  its  volume  exceeded  that  ol  Skylab  1  itsell.  During  its  30-month 
orbital  life,  which  ended  on  11  January  1475,  the  rocket  was  the  brightest  ot 
artificial  satellites  and  on  an  overhead  trunsit  was  the  most  brilliant  object  in 
the  night  sky,  in  the  absence  of  Venus  and  the  Moon. 

The  orbit  was  nearly  circular  and  inclined  at  50°  to  the  equator.  Initially, 
the  uveruge  height  was  400  km  and  the  orbital  period  43.5  minutes.  The  orbit 
seemed  worthy  of  analysis  because  numerous  observations  were  available,  and 
because  the  orbit  passed  very  slowly  through  the  condition  of  31:3  resonance, 
when  the  satellite  makes  31  revolutions  while  the  Earth  spins  twice  (relative  to 
the  orbital  plane).  Analysis  of  the  orbital  perturbations  at  ih;  sonunee 
offered  the  chance  to  obtuin  values  of  lumped  31st-ordor  harmonics  more  accurate 
than  the  only  others  previously  obtained^,  from  l’roton  4,  which  was  not  in  such  a 
favourable  orbiL.  In  addition,  the  variation  in  inclination  away  from  resonance 
was  expected  to  give  values  for  the  rotation  rate  uf  the  upper  atmosphere  over  a 
wide  range  of  heights;  and  there  was  the  possibility  of  testing  the  theory  for 
the  contraction  of  near- circular  orbits  in  an  atmosphere  with  strong  day-to-night 
variation  in  density. 

2  THE  OBSERVATIONS 

During  its  first  15  days  in  orbit,  1473-3711  passed  most  northern-hemisphere 
observing  stations  in  the  daylight  hours,  and  (apart  from  a  sighting  by  .1.  Hewitt 
and  the  author  on  15  May  from  Lagonissi  in  tlreece)  the  only  observations  avail¬ 
able  were  two  from  Sydney,  Australia,  and  observations  on  8  transits  from  the 
radar  at  RSRE  Malvern  -  not  enough  to  allow  the  determination  of  unbiased  orbits. 

From  June  1473  onwards,  daily  US  Navy  observations  from  the  Navspasur 
system  were  available*,  and  orbits  could  have  been  determined  at  any  chosen  time. 
Continuous  otbit  determination  was  not  necessary  for  1473-3711,  however,  except  at 
the  time  of  the  31:2  resonance  and  at  the  end  of  the  lile,  so  the  policy  adopted 
was  to  determine  orbits  whenever  there  were  appreciable  numbers  ol  other  observa¬ 
tions  in  addition  to  the  US  Havv  observations.  About  1500  ot  the  US  Navy 
observations  were  used,  while  the  'oliiei  obsci \ j, i oos 1  tell  into  live  categories, 
as  follows.  _ _ _  _  __  _  _ _ _ 

*  Navspasur  observations  were  being  made  in  May,  ol  course,  but  were  not 
requested  until  June. 
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First,  and  most  accurate,  were  28  transits  recorded  by  the  Hewitt  cameras 
at  Malvern  (27  transits)  and  Edinburgh  (1  trunsit),  with  accuracies  between  1 
and  3  seconds  of  arc.  A  maximum  of  5  observations  per  transit  was  used,  although 
on  some  transits  only  one  or  two  observations  were  available. 

Second  in  order  of  accuracy  was  the  kinetheodolite  at  the  South  African 
Astronomical  Observatory.  These  observations,  of  accuracy  between  30  and  100 
seconds  of  arc,  were  available  on  69  transits,  usually  with  2  observations  per 
transit,  and  had  a  profound  influence  on  the  values  of  the  orbital  parameters, 
because  of  the  accuracy  of  the  observations  and  the  sparsity  of  other  southern- 
hemisphere  data. 

The  third  and  largest  group  of  'other  observations',  over  1400,  came  from 
volunteer  visual  observers  reporting  to  the  Appleton  Laboratory  at  Slough  or  the 
Moonwatch  Division  of  the  Smithsonian  Astrophysical  Observatory.  For  the 
majority  of  these  observations  the  accuracies  were  between  I  and  4  minutes  of  arc 
cross-track. 

The  fourth  group  of  observations  was  from  the  radar  tracker  at  RSRE  Malvern. 
Observations  were  available  on  82  transits  and  a  muximuw  of  3  observations  per 
trunsit  was  used,  although  sometimes  only  one  was  available.  The  average 
accuracy  was  about  1  km  in  range  and  3  minutes  of  arc  in  direction. 

The  fifth  group  consisted  of  more  than  1000  observations  in  the  last  7  days 
of  the  satellite's  life,  provided  by  the  assigned  and  contributing  sensors  of  the 
North  American  Air  Defense  Command  (NORAD)  Space  Detection  and  Trucking  System 
(SPADATS).  With  these  observations,  and  the  UP  Navy  and  Malvern  radar  observa¬ 
tions,  it  was  possible  to  determine  0  excellent  orbits  tit  daily  intervals  at  the 
end  of  the  life. 

3  THE  ORBITS 

3 .  1  General 

3  4 

The  orbits  were  determined  using  the  l’ROI’6  orbit  refinement  program  ’  ,  in 
which  each  observation  is  given  an  a  priori  error,  and  the  program  proceeds  by  an 

iterative  process  to  reduce  the  measure  of  lit,  >  ,  where  t ^  is  defined  as  the 

sum  of  the  squares  of  the  weighted  residuals,  divided  by  the  number  ol  degrees  of 
tiecdoiu.  The  weighted  residual  is  the  residual  in  light  ascension  or  declination, 

divided  by  the  a  priori  error.  If  the  orbit  fits  the  observations  well,  the 

value  of  c  should  be  somewhat  less  than  1,  because  the  ii  prim  errors  are 
normally  taken  as  round-number  values  somewhat  higher  than  the  actual  errors.  In 
the  orbit  refinement  process,  any  observations  with  weighted  residuals  greater 
than  Ji  arc  rejected. 
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Thu  orbit  of  I973-27B  proved  unexpectedly  difficult  to  determine,  chiefly 
because  an  excessive  number  of  observations  were  rejected.  Two  reasons  for  this 
difficulty  are  now  apparent.  First,  the  brightness  of  1973-2711  enticed  into 
action  a  large  number  of  new  visual  observers,  whose  accuracy  did  not  match  their 
enthusiasm;  in  retrospect  it  would  have  been  more  efficient  not  to  have  used  any 
of  their  observations.  (The  experienced  visual  observers  gave  realistic  esti¬ 
mates  of  their  accuracy,  as  usual.)  Second,  nearly  all  the  observations  were 
less  accurate  than  usual,  because  of  the  satellite's  large  size,  great  brightness 
and  high  apparent  angular  rate  of  travel  (see  section  4),  Over  all,  about  252  of 
the  observations  were  rejected. 

3 , 2  Accuracy  of  the  orbital  parameters 

Orbits  were  determined  at  62  epochs  between  17  June  1973  and  11  January 
1975.  In  the  first  year,  before  the  31:2  resonance  began  to  take  effect,  orbits 
were  determined  whenever  there  were  substantial  numbers  of  other  observations  to 
support  the  US  Navy  observations:  orbits  were  calculated  at  27  epochs  during 
these  12  months,  and  US  Navy  observations  at  elevations  less  than  20°  were  not 
used.  During  the  months  when  resonance  effects  were  appreciable  (1  July  to 
10  October  1974),  as  many  orbits  as  possible  (17)  were  determined,  and  the  low 
elevation  US  Navy  observations  were  included.  Between  16  October  1974  and  decay, 
the  procedure  reverted  to  the  less  intensive  mode,  and  18  orbits  were  determined, 
including  the  6  daily  orbits  in  the  last  6  days.  The  values  of  the  orbital  para¬ 
meters  are  listed  in  Table  1  (pages  43-45). 

The  accuracy  of  the  orbits  varied  a  great  deal,  the  presence  or  absence  of 
Hewitt  camera  or  kinetheodolite  observations  being  a  decisive  factor.  On  16  of 
the  18  orbits  with  Hewitt  camera  observations,  the  standard  deviation  in 
inclination  i  was  between  0.0001°  and  0.0007°,  with  a  mean  of  0.0004° 

(equivalent  to  50  m  in  position).  On  15  of  these  18  orbits,  the  standard 

_6 

deviation  in  eccentricity  e  ranged  between  4  and  15  *  10  ,  with  an  average  of 

10  x  10  equivalent  to  70  m  in  perigee  distance;  the  standard  deviation  in 
right  ascension  of  the  node  if  was  0.002°  or  less  on  all  18  orbits. 

On  40  of  the  44  orbits  without  Hewitt  camera  observations,  the  standard 
deviation  in  inclination  ranged  between  0.0005°  and  0.0016°,  with  an  average  of 
0.0010°,  equivalent  to  120  m  in  distance.  (The  remaining  four  orbits,  for  which 
observations  were  deficient  in  either  number  or  variety,  had  accuracies  between 
0.0018°  and  0.0025°.)  On  40  of  the  44  orbits,  the  standard  deviation  in 
eccentricity  rangeu  from  5  to  21  x  10  with  a  mean  of  12  *  10  equivalent  to 
80  m  in  perigee  height;  and  the  average  standard  deviation  in  11  was  less  than 
0.002°. 
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Thus  the  presence  of  Hewitt  camera  observations  improves  ihe  sd  in 
inclination  by  a  factor  of  2.5,  but  in  eccentricity  by  a  factor  of  only  1.2:  tor 
accurate  values  of  eccentricity,  a  wide  geographical  distribution  of  observations 
is  more  important  than  very  accurate  observations  from  one  site,  and  the  Cape 
kinethcodolite,  because  of  its  geographical  advantage,  was  often  just  as 
influential  as  the  Hewitt  camera  in  determining  e  . 

Within  PkUP,  the  mean  anomaly  M  is  expressed  as  a  polynomial  in  time  t  , 
of  the  form 

M  -  M0  ♦  M  t  t  >!.,t2  +  ...  M  f’  ,  ( 1 ) 

where  t  is  the  time  in  days  after  epoch,  Mq  is  the  value  of  M  at  epoch,  and 
M,  ...  are  coefficients  which  are  included  if  they  are  needed.  With  ll>  of 

the  orbits,  only  M^t  Mj  and  M.,  were  required:  this  implies  that  the  drag 
remained  nearly  constant  during  the  time  interval  covered  by  the  observations. 

On  the  remaining  orbits  larger  numbers  of  M-coef f icients  were  needed,  as  a  result 
of  variations  in  density.  The  density  at  heights  near  240  km  during  1973  and 
1°74  has  been  determined  by  Walker^  at  average  intervals  of  about  2  days: 
irregular  day-tr — day  variations  of  up  to  about  20”  often  occur,  and  the  amplitude 
of  such  variations  would  be  larger  at  the  heights  being  sampled  by  1 ^ 7 3—27 H 
during  its  first  year  in  orbit  (350-400  km).  The  density  at  these  heights  at 
weekly  intervals  during  1974  has  been  determined  by  Brookes  and  Moo re ^  from 
analysis  of  the  orbit  of  i973-27A. 

The  use  of  numerous  M-cocff icients  is  not  to  he  recommended,  because  of 
their  possible  correlations  with  other  parameters.  Here  on  extra  M-coef T ioient 
has  been  used  only  if:  (a)  the  value  of  the  measure  of  fit  r  is  reduced  by  at 
least  20”;  (b)  the  standard  deviations  of  the  main  orbit  parameters  are  reduced; 
and  (c)  the  value  of  the  new  M-coef f ioient  is  at  least  5  times  its  standard 
deviation.  The  ingredients  of  a  good  orbit  -  exemplified  by  No.  S  -  are 

(i)  having  accurate  observations,  well  distributed  geographically,  and 

(ii)  fitting  with  Mq,  and  only,  preferably  with  a  low  value  of  M.,  . 

The  values  of  M  and  M.,  are  of  course  very  we  11 -determined  on  all  the 
*  *■  ^ 

orbits.  Mj  being  accurate  to  1  part  in  10  (except  on  the  last  orbit)  and  M, 
to  1  part  in  100  (except  on  orbits  4,  26  and  28). 

For  most  of  the  orbits,  changes  in  the  assumed  value  for  the  density  scale 
height  H  had  a  negligible  effect,  and  H  was  taken  as  50  km.  The  effect  of  H 
was  appreciable  only  on  the  last  four  orbits,  when  the  perigee  heights  were 


7 


.'01,  193,  180  and  158  km  respectively,  and  values  of  H  were  taken  as  HI,  Jl>,  J(< 
and  .’0  km  respectively,  in  accordance  with  the  COSl’AR  International  Reference 
Atmosphere  1972  (Ref  7). 

Table  1  shows  the  wide  variation  in  the  number  of  observations  per  orbit, 
between  20  and  98,  with  a  mean  of  45.  The  time  intervals  covered  by  the  observa¬ 
tions  range  between  1  day  (for  the  last  6  orbits)  and  9  days,  with  a  mean  of  5 
days.  The  values  of  e  (except  on  orbit  57)  were  between  0.40  and  0.88,  which 
is  very  satisfactory. 

The  last  six  (daily)  orbits  were  of  excellent  accuracy  and  consistency.  The 
satellite  burnt  up  near  the  Azores  at  07.05  UT  on  11  January  1^75,  so  the  epoch 
of  the  final  orbit  is  only  7  hours  before  decay.  Two  US  Navy  observations  at 
05.59  UT  were,  not  surprisingly,  rejected  by  TROT,  and  their  sd  was  artificially 
increased  (to  about  10  km)  to  ensure  that  they  were  used  in  the  final  orbit,  to 
strengthen  the  values  of  the  M-coeff icients. 

3.3  Variation  of  the  orbital  parameters 

The  values  of  inclination  on  the  62  orbits  of  Table  1,  with  their  standard 

deviations,  are  plotted  in  Fig  1.  The  general  decrease  due  to  the  effect  of 

atmospheric  rotation  is  visible,  with  a  hint  of  an  increase  at  the  time  of  the 

31:2  resonance;  but  the  values  are  rather  scattered,  and  no  conclusions  can  be 

drawn  until  the  lunisolar,  zonal  harmonic  and  J,  ,  perturbations  are  removed  (see 

*-  *  — 

section  5). 

The  62  values  of  eccentricity  are  plotted  in  Fig  2,  their  standard  devia¬ 
tions  being  too  small  to  indicate.  Fig  2  shows  clearly  the  oscillation  due  to 
odd  zonal  harmonic  perturbations,  which  is  of  amplitude  0.8  x  10  together  with 
the  general  decrease  in  eccentricity  due  to  air  drag.  The  points  have  been 
joined  by  a  broken  curve  as  a  guide  to  the  eye.  Any  perturbation  due  to  the  31:2 
resonance  must  be  very  small  compared  with  the  major  perturbations  due  to  zonal 
harmonics  and  air  drag. 

Fig  3  gives  an  alternative  plotting  of  the  eccentricity,  using  e  and  the 
argument  of  perigee  to  as  polar  coordinates.  In  the  absence  of  drag,  the  locus 
of  the  points  would  be  a  circle  with  its  centre  at  the  point  C  on  the  e  sin  ,o 
axis:  the  action  of  drag  converts  the  circle  into  a  decaying  spiral.  The  day-to¬ 
night  variation  in  air  density  displaces  the  centre  of  the  spiral  from  the 
g 

e  sin  u>  axis  ,  but  this  is  scarcely  perceptible  in  Fig  3,  where  the  spiral  is  of 
classic  regularity  until  the  last  few  days  in  orbit.  (The  variation  of  e  and 
to  in  these  last  days  is  discussed  in  section  8.3.) 
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Kig  4  shows  the  values  of  the  light  ascension  of  the  node,  which  decreases 
smoothly  and  ulmout  linearly,  at  a  rate  of  S ,  2  deg/dav  initially,  increasing  to 
‘'.d  deg/day  on  the  final  orbit. 

Fig  3  gives  the  values  of  the  orbital  decay  rate  M,  .  As  time  gasses,  the 
satellite  gradually  sinks  lower  into  the  atmosphere,  and  there  is  a  general 
increase  in  the  decay  rate.  The  values  also  reflect  the  semi-annual  variation  in 
air  density'*,  with  minima  of  M.,  in  August  1973  and  January  1974,  separated  by 
maxima  in  October  1973  and  March  1974.  The  minimum  in  July  1974  is  just  visible, 
but  the  massive  increase  in  September-Oe tuber  1974  is  caused  mainly  by  the 
decreasing  height  of  the  satellite. 

4  ANALYSIS  OF  THE  OBSERVATIONS 

The  residuals  of  the  observations  utilized  in  orbits  l-3ti  have  been  listed 
with  the  aid  of  the  ORES  program  ,  and  have  been  sent  to  the  observers.  Table  2 
gives  the  nns  residuals  for  the  stations  (in  nine  different  countries)  having 
five  or  more  observations  accepted.  The  angular  residuals  tor  Stution  29  are 
geocentric,  and  need  to  be  multiplied  by  a  factor  of  about  3  to  make  them 
comparable  with  the  other  (topoeontrie)  values.  The  residuals  for  observations 
on  the  lust  six  orbits  are  not  given. 

The  residuals  in  Table  2  are  generally  rather  larger  than  in  other  recent 
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orbit  determinations,  using  similar  observations  .  One  likely  reason  is  that 

1973-27B  often  moved  rapidly  (at  1  dog/s  or  more)  relative  to  the  star  back¬ 
ground;  such  rapid  motion  increases  the  angular  error  resulting  from  timing 
errors.  The  Hewitt  camera  at  Malvern,  which  has  vury  accurate  timing,  is 
unaffected  by  the  rapid  motion,  and  the  total  nns  residual  for  the  camera  is 
only  2  seconds  of  arc,  as  usual.  But  for  the  visual  observers  the  timing  errors 
when  observing  1973-2711  (about  100  ms)  were  more  important  than  angular  errors, 
and  the  along-track  residuals  usually  exceeded  the  cross-track  residuals.  This 
tendency  is  visible  in  Table  2  for  tbo  British  observers,  who  often  made  their 
observations  of  1973-2711  when  it  was  at  high  elevation  in  the  south  and  was 
travelling  from  west  to  east:  in  these  circumstances  the  along-track  errors  are 
mainly  in  right  ascension,  and  the  tight  ascension  residuals  are  on  average  about 
twice  those  in  declination  for  the  seven  British  observers  with  the  most  observa¬ 
tions  accepted  -  Stations  2420,  2414  ,  2421  ,  22HS,  24I9(  2430  and  2428.  The 
declination  residuals  therefore  give  a  better  idea  of  the  angular  accuracy  ot  the 
observers,  and,  since  the  mis  residuals  tend  to  be  unduly  influenced  by  one  poor 
observation,  it  is  fair  'o  say  that  in  good  conditions  the  visual  observers 
should  he  capable  of  an  angular  accuracy  about  half  the  nns  declination  residuals 


him . . 


awnr-ai; 
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in  Table  2.  (The  relation  between  nns  residuals 
conditions  has  been  discussed  in  section  5  of  Kef 


Table  2 

Residuals  of  the  observations  on  orbit 


with  5  or  more  observations 


Station 


Number  of 
observations 
accepted 


1  US  Navy 

2  US  Navy 

3  US  Navy 

4  US  Navy 

5  US  Navy 

6  US  Navy 

29*  US  Navy 

86  4  2597  Rodcntal 

161  Castlcrock 

414  Capetown 

514  Thai  Ion 

7 l 9  Rodewisch 

759  Cluj 

265  Farnham 

303  Malvern  camera 

04  Malvern  radar 

14  Bournemouth 

19  Tremadoc 

20  Willowbrae 

21  Malvern  4 

28  Fr imley 

30  Stevenage  4 

28  Aldershot 

39  Dymchurch 

73  Genoa  1 

77  Cape  kinetheodolite 

96  Akrotiri 

26  Groningen 

30  Denckamp 

56  Apeldoorn 

59  Ache  1  2 

65  Groningen  3 

17  Sacramento  I 

Others 


Geocentric 


and  the  average  in  best 
13.) 


s  1-56,  for  stations 
accepted 


Rms  residuals 

Range 

Minutes  of 

arc 

km 

RA 

Dec 

Total 

2.2 

2.3 

3.2 

2.5 

3.0 

3.9 

1  .7 

1  .6 

2.3 

1  .6 

1.9 

2.5 

2.  1 

2.  1 

3.0 

2.2 

2.3 

3.2 

0.6 

0.4 

0.4 

5.  1 

3.0 

5.9 

1  .8 

1 .  1 

2.1 

3.3 

4.  1 

5.3 

3.9 

4.0 

5.6 

5.0 

3.4 

6.0 

1  .8 

4.7 

5.0 

4.0 

2.  1 

4.6 

0,03 

0.02 

0.04 

1  .3 

4.4 

2.8 

5.7 

3.0 

6.4 

4.8 

2.6 

5.4 

4.6 

1  .  5 

4.8 

3.6 

2.5 

4.4 

6.0 

3.3 

6.8 

4.0 

2.0 

4.5 

2.6 

3.0 

4.0 

3.5 

1.8 

3.9 

8.6 

16.1 

18.3 

1  .2 

1.3 

1  ,1 

4.2 

4 . 6 

6.2 

2  2 

1  .5 

2.7 

4.8 

4.7 

6.7 

5.7 

3.5 

6.7 

4.3 

6.7 

8.0 

5.2 

5.3 

7.4 

5.3 

4.8 
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It  is  also  probable  that  the  accuracy  of  the  visual  observers  was  adversely 
affected  by  the  extreme  brightness  of  the  satellite,  which  was  sometimes  of 
magnitude  -2.  Reference  stars  were  ol Leu  losL  in  the  glare  as  the  satellite 
passed  them,  a  situation  scarcely  conducive  to  good  accuracy.  This  difficulty, 
aggravated  by  the  high  angular  velocity  of  the  satellite  (often  more  than  1°  per 
second),  proved  too  much  for  some  of  the  inexperienced  visual  observers,  and  only 
743  of  the  original  1400  visual  observations  were  accepted  in  the  final  orbits. 

Of  these,  18%  were  from  Willowbrae  (R,D.  Eberst),  15%  from  Bournemouth 

(D.J.  Hopkins),  10%  from  Malvern  4  (D.M.  Brierley),  6%  from  Cape  Town 

(W.P,  Hirst),  5%  from  Farnham  (D.G.  King-Hele)  and  4%  from  Tremadoc  (D.C.  Mason). 

Thus  58%  of  the  accepted  visual  observations  came  from  only  6  observers  -  a 
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pattern  familiar  m  many  previous  orbit  determinations 

The  US  Navy  and  Malvern  radar  observations,  for  which  timing  is  by 
electronic  methods,  also  have  rather  larger  residuals  than  usual:  probably  the 
errors  were  increased  because  the  satellite  was  so  large  that  it  produced 
multiple  radar  reflections. 

The  kinetheodolite  observations  were  also  less  accurate  than  usual.  There 
are  two  possible  reasons:  first,  the  unusually  rapid  angular  tracking  rates 
needed  on  transits  at  high  elevations  might  have  had  a  bad  effect;  second,  and 
rather  paradoxically,  the  clear  skies  of  South  Africa  and  the  brightness  of  the 
satellite  combined  to  allow  an  exceptional  number  of  observations  at  very  low 
elevation,  often  less  than  10°,  for  which  errors  in  the  atmospheric  refraction 
correction  would  be  larger  (see  section  5) . 

5  REVISION  OF  THE  ORBITS 

When  lunisolar,  zonal  harmonic  and  J?  ?  perturbations  have  been  removed,  the 

*■9 

values  of  inclination  (outside  the  31:2  resonance  region)  should  lie  on  a  smooth 

curve,  slowly  decreasing  under  the  influence  of  atmospheric  rotation.  (The 

14 

effects  of  other  perturbations,  such  as  Earth  and  ocean  tides  ,  are  expected  to 
be  less  than  the  standard  deviations  of  the  values,  and  are  ignored.) 

With  the  orbits  originally  computed,  however,  the  values  of  inclination, 
after  removal  of  perturbations,  did  not  lie  on  a  smooth  curve,  but  were  quite 
scattered.  The  theoretical  curves  for  various  atmospheric  rotation  rates  A 
were  calculated,  and  a  value  A  =  1.0  between  launch  and  resonance  seemed  to  fit 
the  scattered  values  best:  after  some  tr ial-and-error  experiments  with  orbits 
1-7,  a  mean  curve  was  established,  with  A  =  1.0  .  On  the  assumption  that:  the 
mean  curve  was  more  likely  to  give  the  correct  value  of  i  than  the  scattered 
orbits,  any  orbit  between  No. 8  and  No. 27  that  did  not  fit  the  curve  was  re-run 
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with  the  inclination  fixed  at  the  value  indicated  by  the  curve.  Any  observations 
which  could  not  tolerate  this  inclination  were  then  omitted  -  with  the  exception 
of  the  Cape  kinetheodol i te  observations,  for  which  the  a  priori  error  was 
doubled  if  the  observation  was  rejected  or  had  a  large  residual.  The  orbit  was 
then  again  re-run,  with  the  inclination  free.  On  most  occasions,  the  revised 
orbit  thus  obtained  had  both  a  lower  value  of  t.  and  lower  standard  deviations 
than  the  original  orbit,  as  well  as  an  inclination  closer  to  the  mean  curve:  this 
is  a  strong  indication  that  the  non-conforming  observations  were  in  error.  (If 
c  was  not  reduced,  the  new  orbit  was  not  accepted.)  A  similar  procedure  was 
applied  af ter  .  resonance .  It  is  of  course  these  revised  orbits  that  are  given  in 
Table  1  and  Fig  1,  and,  when  the  perturbations  are  removed,  the  values  of 
inclination  decrease  quite  smoothly,  as  Fig  6  shows,  except  near  resonance. 

It  may  seem  disquieting  that  such  refinement  is  needed:  the  need  arises 
because  of  the  vulnerability  of  least-squares  fitting  to  the  influence  of 
observations  which  are  in  error  by  two  or  three  times  their  pre-assigned  accuracy. 
Omitting  such  an  observation  can  change  the  value  of  i  by  several  standard 
deviations.  'Wrong'  observations  are  often  very  difficult  to  identify,  however  - 
especially  if  isolated  geographically  -  when  fitting  a  changing  orbit  to  the 
motion  of  a  satellite  decaying  under  the  action  of  variable  air  drag.  To  reduce 
the  rejection  level  to  less  than  3c  would  be  to  risk  rejecting  good  observa¬ 
tions,  and  a  rejection  level  differing  between  one  orbit  and  another  would  be 
difficult  to  justify.  So  the  technique  of  intr  .-ucing  ir.deper,'Je 't  information 
about  the  likely  value  of  inclination  seems  .to'  tnly  legitimate  and  logical,  but 
quite  essential,  if  the  dilemma  is  to  be  resolved*. 

For  1973-27B,  the  problem  can  be  seen  in  retrospect  as  largely  stemming 
from  overweighting  of  the  Cape  kinetheodoli te  observations.  These  have  a  greater 
influence  than  any  other  observations,  because  they  are  often  both  geographically 
isolated  and  more  accurate  than  any  others  (in  the  absence  of  Hewitt  camera 
observations).  In  the  original  orbits  the  accuracy  of  the  kinetheodoli te 
observations  was  taken  as  0.01°  whenever  this  was  the  figure  specified  by  the 
observers.  It  would  have  been  wiser  to  have  used  0.02°,  for  several  reasons. 

*  Though  the  least-squares  method  has  its  virtues,  orbit  determination  tends  to 
expose  its  weaknesses.  Godwin  said  of  Parliament,  To  decide  truth  by  the 
casting  up  of  numbers  is  an  intolerable  insult  to  reason  .  To  decide  that 
'truth'  is  revealed  by  the  exact  location  of  a  shallow  minimum  in  a  nearly 
flat  multidimensional  surface  may  also  seem  an  insult,  when  the  removal  of 
only  one  observation  can  so  alter  the  surface  that  the  minimum  may  be 
displaced  by  an  amount  much  greater  than  the  formal  sd.  It  seems  unlikely 
that  other  methods  can  do  much  better,  however,  because  there  is  no  sure  way 
of  defining  a  'bad'  observation. 


12 


First,  although  0.01°  is  a  fair  ' round-number '  accuracy  estimate,  the  average 
error  of  those  observations  is  probably  rather  larger,  perhaps  about  0.013°. 
Second,  all  the  other  observations  ure  given  pessimistic  a  pi'iot'i  accuracies  - 
that  is  why  the  values  of  l  are  usually  between  0.5  and  0.8  -  and  it  is  bud 
practice  to  give  heavier  relutive  weighting  to  isolated  observations.  Third, 
there  are  some  doubts  about  tie'  accuracy  of  the  station  coordinates  for  the 
kine theodolite,  which  would  argue  for  an  error  larger  than  0.01°  for  high- 
elevation  observations,  Fourth,  us  mentioned  in  section  4,  many  of  the  kine- 
thoodolite  observations  of  1 9 7 3 — 2 7 B  were  made  at  very  low  elevutions,  where  the 
larger  errors  in  the  refraction  correction  might  vitiate  the  assumption  of  an 
accuracy  of  0.01°.  For  these  reasons,  it  is  roeonmended  that  an  uccuracy  of 
0.02°  should  be  used  tor  the  Cape  kinetheodolite  observations  in  future  orbit 
determinations:  this  will  ensure  that  tuturo  observations  will  not  be  lost 
through  rejection  (a  loss  thaL  can  be  ill  afforded);  and  fewer  orbits  will  he 
disturbed  by  the  minority  of  observations  that  would  otherwise  have  been 
overweighted. 

b  THE  31:2  RESONANCE 

b. 1  General 

The  orbit  of  1973-278  experienced  the  effects  of  31:2  resonunce  between 
July  and  October  1974,  Exact  31  :2  resonance  occurred  on  23  Augvist  1974 
(MJ1)  42282),  The  ruw  values  of  inclination  in  Fig  l  give  a  hint  of  an  increase 
in  inclination  ut  the  time  ot  the  resonance,  and  this  emerges  more  clearly  when 
the  perturbations  are  removed,  as  in  Fig  b .  Any  variations  in  eccentricity  due 
to  resonance  are  masked  in  Fig  2  by  the  larger  variations  due  to  odd  zonal 
harmonics  and  air  drag.  The  removal  of  these  effects  is  difficult,  and  is 
discussed  in  section  b.4, 

b . 2  Theory 

The  longi tude-depemlent  part  of  the  Earth's  gravitational  potential  at  an 
exterior  point  (r,  0,  \)  may  be  written  in  normalized  form* J  ns 


,  s in  iiioN,  ,  (.2) 

tin  )  fin 

wh»re  r  is  the  distance  from  the  Earth's  centre,  e  is  co-latitude,  \  is 
longitude  (.positive  to  the  east),  r  is  the  gravitational  constant  lor  the  Faith 

ni 

(  i'WiOti  km  /s  ),  It  is  the  Fai  th's  equatorial  radius  ((>378.  1  km),  1’  (cos  '*)  is 


l-1''(cos  n ) 


{c 

\  Vm 


cos  m  V  +  S 


m^  1 
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the  associated  Legendre  function  of  order  m  and  degree  V  ,  and  C  and  Sfin 
are  the  normalized  tessera!  harmonic  coefficients.  The  normalizing  factor  N 


is  given  by 


2(21  +■  1  )  (t  -  m)  1 
(&  *  m)  ! 


When  a  satellite  orbit  experiences  fl:«  resonance  -  that  is,  when  the 
satellite  makes  ft  revolutions  while  the  Earth  makes  a  revolutions  relative  to 
the  satellite’s  orbital  plane  -  the  rates  of  change  of  inclination  and  eccen¬ 
tricity  depend  on  the  resonance  angle  ,  defined  by 

<!■  -  u(u)  +  M)  +  |i(S  -  v)  ,  (4) 

where  v  is  the  sidereal  angle.  The  perturbation  produced  in  i  and  e  by  a 
relevant  pair  of  harmonic  coefficients,  and  ,  may  be  written**’'*'7 


di  =  n (  1  -  e  ) 
dt  sin  i 


(l)  hm,.Ct|,„(k  cos  i-»).*|jL"n'*l(cta-isJcxp|j(, «-,«>)  ,  (5) 


o  I 

\a)  tmp  2.pq} 


|^[J  l  "m+  1  ^Sra  '  J  eXl>  1 J  (y*  ' 


....(6) 


Here  !•’„  is  Allan's  normalized  inclination  function  ,  (1,  is  a  function  of 

£mp  _  fpq  17 

eccentricity  e  for  which  explicit  forms  have  been  derived  by  Gooding  , 
denotes  'real  part  of'  and  j  ■  /-l  .  The  indices  y,  p,  k.  and  p  are  integers, 

with  y  taking  the  values  1,  2,  3,  and  q  the  values  0,  ±1,  ±2 .  The 

equations  linking  f,m,k  and  p  are*  :  m  =  yfi;  k  =  ya  -  q;  2p  ■  1  -  k  . 

For  a  particular  H:n  resonance,  the  relevant  value  of  m  in  the  , 

S,  coefficients  is  decided  by  the  choice  of  y  ,  since  m  *»  y|l  .  There  is  a 

series  of  possible  values  of  f  -  any  values  such  that  t  S  m  and  (t  -  k)  is 

even.  The  successive  C„  ,  S,  coefficients  th;it  arise  (for  given  >  and  q) 

«m  vm 

may  usefully  be  combined  in  'lumped'  form  and  written  as 


l.k  V"*  ll’k- 

,  =  ’ 


Y'  fi.>- 
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i 

whore  i  increases  in  steps  of  2  from  its  minimum  permissible  value  tf  ;  tile 

ok....  ^ 

Q,  are  functions  of  inclination  that  can  be  taken  as  constant  for  a  parti- 

o  k 

culur  satellite;  and  Q -  1  when  t  -  <  . 

Here,  with  |5  -  31  and  u  -  2  ,  the  relevant  harmonics  arc  of  order 
m  -  3 1 Y  ,  where  i  -  1,  2,  3,  ...  .  The  harmonics  of  order  02,  93,  ...,  are 
unlikely  to  have  much  influence,  and,  all  hough  their  effects  need  to  be  tested, 
only  the  y  -  1  terms  (that  is,  the  effects  of  harmonics  of  order  m  »  31)  will 
be  given  explicitly. 

Also,  since  k  -  2  -  q  if  \  -  1  ,  the  values  of  (C-  k)  are  even  if 
(if  +  q)  is  oven:  so  the  relevant  values  of  t  are  those  for  which  (C  +  q)  is 
even  and  C.  Z-  31  . 


As  pointed  out  by  Walker*^,  the  q  *  0,  4  1 ,  +2,  ... 
expression  (31  for  di/dt  have  ('.  coefficients  ol  order 


,  terms  in  the 
1  >  s  ( v’ l' )  "  >  •  • 


I  7 

since  the  (> ,  are 
C  pq 

so  j  V o  «•  0,1  for  t  • 
variation  of  i  .  The 
because  in  equal  ion  (h) 


({<*•) 

of  order  - - .  lor  1973-2715,  e  is  ol  order  0.03, 

(hi)! 

bO  ,  and  only  the  q  -  0  terms  are  likely  to  affect  the 

situation  is  different  with  the  perturbation  in  e  , 

the  term  -  II,  iq  -  4  (k  *  lOe^l  is  of  order  e  lor 
e  5  p<|  (  1  *  / 


1  1  •— 

q  ■  0  ;  of  ordcs-^yV  for  q  -  M  ;  of  order  ji'V  for  q  *>  4  2  ;  of  order 
3  2 

<  e  for  q  »  *3  ;  ...  .  With  e  0.005  ,  only  the  q  "  M  terms  are  likely 
to  affect  the  variation  ol  eccentricity. 


So,  when  fitting  theoretical  curves  to  the  variations  of  i  and  e  at 
resonance  lor  1973-2715,  it  is  probable  that  the  only  terms  needed  will  he  those 
with  (y,  q)  **  (1,  0)  for  inclination,  and  (> ,  q)  -  (1,  1)  and  (1,  -I)  for 
eccentricity.  Explicit  forms  for  di/dt  and  de/dt  in  terms  of  these  three 
(Y,  q)  pairs  have  been  given  in  Kef  2.  Here  it  is  only  necessary  to  quote 
explicitly  the  (> ,  q)  *  (1,  0)  terms  for  di/dt  , 


di 

ilt 


COM  >)l'32i3,  t  15jS:1 


31 


sin  T  +  0,^  cos  >1 


*  terms  in 


cos 


(l't  -  q»') 


(K) 
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and  the  (y,  q)  “  (I,  *1)  terms  for  de/dt  , 


de 

dt 


fR\31  ,7S  Ir'*1 

n\a  J  '  17I,31 .31.1 5)  31 


_ltl 

CQ1  sin($  -  d)  -  cos  (l1  -  ui) 


(--1,3  _-l,3 

+  13F3]  3j  |4jC3i  sin(4>  +  w)  -  S ^ j  cosO  +  id) 


+  terms  in  e^l  '  (q  -  ye*1)  “  qw) 


2.  cos  1 


(9) 


In  equations  (8)  and  (9),  the  value  of  >1'  is,  from  equation  (4),  g i von  by 


2 (u<  +  M)  +  31  (12  -  v) 


(10) 


The  values  of  $  and  <1>  for  1973-27B  between  July  and  October  1974  are  plotted 
in  Fig  7, 

The  lumped  C  and  S  coefficients  in  equations  (8)  and  (9)  may  be 

written  in  terms  of  the  F.  und  R/a  as  follows: 
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and  similarly  for  S  ,  on  replacing  C  by  S  throughout.  The  numerical  values 
of  the  lumped  coefficients  for  1973-27)1,  at  inclination  50.04°,  are 


C31  "  ^32,31  7 * A^34 , 3 1  +  /7'6S6,3l  "  4C38, 3 1  +  89,0^40,31  /1,l842,31 


<r 

o 


+  7 "  "^44 , 3 1  +  A5,2^46,31  "  33'iC48,31  ~  l7‘IC50,31  +  3I*^52,; 


+  3.0C  .  ^  -  24.60  ,,  +  2.50  u  +  18. 40  -  .... 

54,31  56,31  58,31  60,31 


31 


....U4) 
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and  similarly  for  S  .  liquations  (14)  —  < 16)  show  that,  if  good  values  for  the 
lumped  coefficients  can  be  found,  they  will  provide  a  stringent  test  of  the 
accuracy  of  future  geopotential  models  with  coefficients  to  degree  60  or  more. 

6,3  Analysis  of  inclination 

The  analysis  begins  50  days  before  exact  resonance,  at  MJD  42232  when 
I  =  -26.4  deg/day,  and  ends  48  days  beyond  resonance,  at  MJD  42330,  when 
!  =  42.9  deg/day.  As  Fig  7  shows,  $  increases  more  rapidly  after  resonance 
than  before,  because  of  an  increase  in  the  orbital  decay  rate,  resulting  from  the 
satellite's  decreasing  height  and  the  semi-annual  increase  in  air  density  between 
August  and  October. 

The  values  of  inclination  were  cleared  of  zonal  -harmonic  and  lunisolar 

1  8 

perturbations  by  using  the  PROD  computer  program  with  1-day  integration  step; 

the  perturbations  due  to  the  J  Lesseral  harmonic,  which  arc  recorded  with 

each  PROP  orbit,  were  also  removed.  The  resulting  values  of  i  in  the  resonance 

1  9 

region  were  fitted  using  the  THROE  computer  program  with  (y,  q)  =  (1,  0).  In 
these  computations  the  value  of  on  orbit  n  ,  at  epoch  t^  ,  was  replaced  by 


because  the  numerical  integration  within  THROE  proceeds  on  the  assumption  that 

M_  remains  constant  between  Lime  t  and  time  L  , .  .  The  use  of  M.  ensures 

l  n  n+ 1 

that  THROE  correctly  models  the  total  change  in  due  to  drag  between  orbit  n 

and  orbit  (n  +  1). 
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TIk>  first  THKOK  run  using  the  17  orbits  with  an  atmospheric  rotation  rate 

A  -  1  .  2  rov/day  cave  encouraging  results:  the  measure  of  fit  t  was  1.03  and  the 
0-0,2  o_0,3 

values  ol  10  and  10  8^  were  0.4  1  O.o  and  -1.0  i  0.4.  However,  several 

ot  the  orbits  fitted  badly,  and  a  process  of  revision  was  started,  similar  to 

that  described  in  section  3.  The  el  feels  ol  changing  the  density  scale  height  11 

and  the  rotation  rate  A  were  also  assessed.  Alteration  of  II  had  little 

effect:  the  value  chosen,  11  «  45  km,  was  that  obtained^  from  CTRA  10?'J  for  the 

appropriate  average  level  of  solar  activity,  and  average  local  time,  at  a  height 

30 

of  350  km,  0.411  above  the  mean  perigee  height  .  Tor  the  rotation  rate  A  , 

3  1 

previous  results'"  indicate  a  mean  value  of  1.3  at  a  height  of  350  km,  and  since 
THROE  revealed  a  slightly  (0.410  better  lit  with  A  ■*  1.3  than  with  A  -  1.2  , 
the  higher  value  was  adopted. 

In  the  final  THROE  run,  witli  the  revised  orbits  and  0 ,  q)  -  (1,  0),  the 
values  obtained  for  the  lumped  JlsL-order  harmonics  were: 


.0  c31 


"  0 . 80  f  0.32, 


(,-0.2 

i°  s3! 


-1.21  ±  0.20  . 


(17) 


The  value  of  i  was  0.5b.  The  fitting  ot  the  theoretical  curve  to  the  individ¬ 
ual  values,  alter  removal  of  perturbations  due  to  atmospheric  rotation  and  luni- 
solar  procession,  is  shown  by  the  unbroken  line  in  Tig  8.  The  tit  is  good,  and 
it  is  perhaps  surprising  that  the  standard  deviations  of  the  coefficients  are  not 

smaller:  this  is  mainly  due  to  the  high  correlation  between  them  (0.89).  The 

1  4 

amplitude  ot  the  oscillation  in  i  due  to  Earth  and  ocean  tides  is  expected  to 
be  about  0.0003°  (35  in),  and  since  only  one  ol  the  17  values  has  a  standard 
deviation  as  low  as  this,  it  is  not  necessary  to  make  any  adjustments  tor  the 
tidal  effects. 

The  orbits  were  also  fitted  using  (i  ,  <|)  »  (I,  0),  (1,  1)  and  (1,  -1);  but 
l  scarcely  changed,  and  the  extra  eocflicicnls  were  indeterminate.  A  further 
fitting,  with  (-,  ,  q)  «  (1,  0)  and  (2,  0),  resulted  in  a  4 %  increase  in  i  , 
little  change  in  the  lumped  list -order  harmonics  (though  their  sd  increased  by 
more  than  75%),  and  indeterminate  values  for  the  b2nd-urdor  terms,  which  thorofoio 
appear  to  be  of  no  significance. 

The  standard  deviations  of  the  values  (17)  should  he  reliable:  THROE  works 
on  the  assumption  that  the  orbits  are  independent,  and  this  seems  fair,  because 
there  were  only  two  orbits  lor  which  a  tew  observations  at  the  end  were  re-used 
at  the  beginning  ot  the  next  orbit,  to  make  up  the  numbers.  Errors  due  to  error 
in  the  chosen  value  of  A  are  not  taken  into  aoeount  ,  hut  should  not  be 
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significant  because  a  change  of  0.1  in  A  alters  the  values  of  the 
coefficients  by  less  than  25?  of  their  mi. 

6.4  Analysis  of  eccentricity 

The  eccentricity  of  a  nearly  circular  orbit  like  that  of  1973-27B  suffers 
important  perturbations  from  the  effects  of  (a)  the  odd  zonal  harmonies  in  the 
geopotential  and  (b)  atmospheric  drag:  both  these  perturbations  need  to  be  care¬ 
fully  removed  if  the  effects  of  the  resonance  are  to  be  exposed.  Lunisolar 
perturbations  are  small  and  easily  removed. 

The  perturbations  produced  by  the  effects  of  odd  zonal  harmonics  on  a  drag- 

">  *■> 

free  orbit  of  low  eccentricity  arc  well-known'*";  but  their  removal  from  an  orbit 
suffering  other  unknown  perturbations  (due  to  drag  and  resonance)  raises  many 
difficulties.  The  method  adopted  here  was  as  follows.  The  variation  of 
eccentricity  with  time  on  a  drag-free  orbit,  initially  coinciding  with  that  of 
1973-27B,  was  calculated  by  numerical  integration  at  1-dav  intervals  with  the 

J  g 

aid  of  the  PROD  computer  program  ,  for  a  gravitational  field  with  zonal 
harmonics  up  to  the  20th,  plus  lunisolar  perturbations.  The  variation  in  e  , 
which  mav  be  denoted  by  Ae  ,  gives  the  perturbation  of  eccentricity  on  the 

£rl 

drag-free  (and  resonance-free)  orbit  due  to  zonal  harmonics  and  lunisolar 

perturbations.  In  the  absence  of  drag,  Ae  depends*  on  the  argument  of  perigee 

L  H 

u>  ,  and  it  was  assumed  here  that  the  appropriate  value  of  Ae  to  be  subtracted 

L.  H 

from  the  observed  value  of  e  at  each  epoch  was  the  value  given  by  the  drag- 

free  calculation  not  at  the  epoch,  but  at  the  slightly  different  date  when  u>  on 

the  drag-free  orbit  had  the  same  value  as  on  the  actual  orbit.  This  difference, 

though  never  more  than  1  day,  was  significant  because  the  values  of  w  on  the 

actual  orbits  differed  from  those  on  the  drag-  and  resonance-free  orbits  by  up  to 

3°,  and  a  3°  change  in  «>  changes  Ac.,.,  by  4  o  10  which  is  comparable  with 

1 

the  standard  deviation  of  some  of  the  values.  If  the  actual  values  of  to 
differed  more  widely  from  those  on  the  drag-free  orbits,  this  method  would  be 
open  to  criticism,  but  it  seems  likely  to  be  accurate  enough  here. 

Fig  9a  shows  the  values  of  eccentricity  from  the  PROP  orbits,  e  ,  as 
circles,  and  the  values  of  e  after  removal  of  odd  zonal  harmonic  and 

- - j - 

*  On  a  drag-free  orbit,  Ae  -  e  -  {  c~  -  2eB  sin  to  +  t>~  {  ,  where  e  is  the 

observational  value  of  eccentricity,  and  8  is  a  constant  expressible"  in 
terms  of  the  odd  zonal  harmonic  coefficients,  ,T^,  etc. 

|^8  -  — -"-n —  (-  +•  terms  in  J^,  etc)J.  For  1973-2711,  8  -  802  '  10  It  e 

is  not  too  small,  8  is  the  amplitude  of  the  oscillation  in  e  caused  bv  odd 
zonal  harmonics. 
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lunisolar  perturbations  by  the  method  described  above,  as  triangles. 

(f.  -  e  -  he„„) .  The  approximately  sinusoidal  variation  in  e  ,  indicated  by 
the  broken  line,  has  been  removed  with  apparent  success,  but  there  remains  a 
considerable  variation  in  e  attributable  to  air  drag  and  possibly  resonance*. 

The  values  of  e^  were  now  used  as  input  for  a  THROE  run,  and  (after 
removal  of  air  drag  perturbations  within  THROE)  were  fitted  by  the  theoretical 
equation  (9).  As  in  t  'tting  of  inclination,  H  was  given  its  mean  value  of 
45  km,  since  the  mean  value  of  z  -  ae/H  was  near  0.4,  and  the  scale  height  for 
use  in  the  theory  should  be  evaluated  at  a  height  0,4H  above  perigee  when 
z  ■  0.4  (see  Fig  6.3  of  Ref  23).  This  fitting  of  Cj  by  THROE  proved 
disastrous,  giving  c  -  8,8  :  it  was  clear  that  a  fundamental  error  was  being 
made. 

The  error  was  fairly  obvious  -  a  failure  to  allow  for  the  day-to-night 
variation  in  air  density.  At  heights  near  350  km  with  low  solar  activity,  the 
maximum  air  density  (in  the  afternoon)  exceeds  the  minimum  (at  night)  by  a  factor 
of  about  3.  The  air  drag  model  within  THROE  takes  no  account  of  this  day  to¬ 
night  variation  in  air  density,  which  can  have  a  great  effect  on  a  near-circular 

Q 

orbit  .  To  illustrate  the  effect,  assume  that  e  ■  0.0035  ,  so  that  apogee  is 

45  km  higher  than  perigee  (if  a  s  6700  km);  with  H  -  45  km  ,  the  density  at 

perigee  would  be  2.8  times  greater  than  the  density  at  apogee,  in  the  absence  of 

day-to-night  variation.  Suppose  first  that  perigee  is  at  the  day-time  maximum  of 

density  and  the  apogee  at  the  night-time  minimum:  then  the  perigee  density  would 

exceed  the  apogee  density  by  a  factor  of  about  9.  In  these  circumstances  the 

...  2 

orbit  contracts  like  a  normal  orbit  of  higher  eccentricity,  with  e  decreasing 

linearly  and  proportional  to  the  remaining  lifetime  of  the  satellite.  Thus, 

2 

since  the  -ining  lifetime  is  about  200  days,  e  would  decrease  by  102  m 

20  days.  1  suppose  the  situation  is  reversed,  with  perigee  at  the  night-time 

minimum  and  apogee  at  the  day-time  maximum.  Then  the  air  density  at  apogee  is 

slightly  greater  than  at  perigee,  and  drag  is  likely  to  be  nearly  constant  round 

2 

the  orbit.  In  these  circumstances  e  would  remain  very  nearly  constant,  and 

might  even  increase  slightly,  instead  of  decreasing  by  102  in  20  days. 

*  Odd  zonal  harmonic  perturbations  have  to  be  removed  outside  THROE,  because  the 
effects  of  air  drag  are  treated  within  THROE  by  using  PROP  subroutines  in 
which  the  values  of  eccentricity  are  corrected  for  air  drag  on  the  assumption 
that  they  are  not  seriously  perturbed  by  any  force  other  than  air  drag.  For 
many  orbits,  this  procedure  is  of  adequate  accuracy,  but  for  1973-278  the 
change  in  the  perturbation  in  e  between  successive  orbits  can  be  an 
appreciable  fraction  of  e  itself:  for  example  between  MJD  42309  and  42317, 
where  10^e  =  2.0  ,  the  perturbation  to  10-^e  changes  by  0.4  (see  Fig  9a), 

In  these  circumstances  the  values  of  c  have  to  be  cleared  of  perturbations 
before  being  used  in  THROE. 
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A  genera  1  analytical  theory  for  the  effects  of  the  day-to-night  variation 

24 

in  air  density  on  orbits  with  z  *■  I  in  available  ;  but  when  z  1  the 

y 

results  are  very  complicated,  except  in  certain  special  cases,  low  drag,  or 
constant  sun-perigee  angle,  or  weak,  day-to-night  variation  I'nl'or tunato ly 
I973-27B  does  not  qualify  for  any  of  the  special  categories,  so  an  analytical 
correction  for  the  day-to-night  variation  in  air  density  is  not  available,  and 
the  correction  has  to  be  made  to  eacli  orbit  individually.  The  correction  has 
been  applied  using  a  simplified  theory  given  in  the  Appendix,  which  assumes  that 
the  atmospheric  density  depends  on  the  geocentric  ungulur  distance  from  the  point 
M  of  maximum  density,  taken  to  be  on  the  equator  at  14h  local  time  (.the  time 
indicated  by  I'i.'iM  lt).’Z)  .  The  geocentric  angle  $  between  M  and  perigee  for 

I073-27B  during  the  resonance  phase  is  shown  in  Fig  10,  with  the  local  time  at 
perigee  marked  on  the  curve,  and  M  at  either  1 4h  or  16h  local  time  (the  time 
indicated  in  Jacobin*  s  1977  model*-'’).  In  calculating  the  effects  between  orbit 

n  and  orbit  (n  +  1)  the  mean  value  of  *  between  the  two  orbits  is  taken,  and 

P 

the  correction  to  be  applied  for  the  effect  of  the  day-to-night  variation,  c'  , 
is  found  by  integrating  from  orbit  to  orbit.  The  values  of  e^  for  the  17 

orbits,  shown  in  Fig  9b,  are  far  from  negligible:  the  total  change  is  0.5  x  10  \ 

which  may  be  regarded  as  equivalent  to  3  km  in  perigee  distance;  and  even  the 
2-hour  change  in  the  time  of  assumed  maximum  density  has  an  appreciable  effect. 

With  e^  available,  a  doubly  corrected  value  of  o  ,  given  by 
e^  ■  Cj  +  can  be  formed  and  the  values  of  e^  can  be  fitted  using  THROF. 

The  previous  attempt  to  fit  Uj  hud  been  disastrous  because  the  fitting  was  not 
only  bad  but  ulso  gave  no  hint  as  to  how  improvements  might  be  made:  nearly  all 

the  values  seemed  to  be  wrong  -  as  in  fact  they  were.  When  the  values  of  o  , 

were  fitted  with  equation  (9)  by  using  THROF  there  was  a  great  reduction  in  i  , 
from  8.8  to  4.1,  and  there  was  also  a  clear  indication  that  three  values  (at  MJD 

42232,  42245  and  42303)  were  failing  to  fit.  It  would  have  been  possible  to  omit 

♦ 

these  values,  but  the  numerical  integration  of  >1  would  then  be  over  different 
time-intervals  in  the  fittings  of  i  and  ol  c  :  it  is  better  to  keep  the  values 
of  exactly  the  same  as  on  the  fittings  of  i  and  of  e  ,  by  adopting  the 
alternative  strategy  of  increasing  the  standard  deviations  of  the  ill-fitting 
values  by  an  arbitrary  factor  of  10,  to  destroy  their  influence.  With  this 
alteration,  the  value  of  r  decreased  from  4.1  to  1.72. 

The  effect  of  changing  the  assumed  local  time  of  the  maximum  density  from 
14h  to  !(>h  was  now  tested.  The  value  of  <  increased  from  1.72  to  2.95,  so  it 
was  clear  that  maximum  density  at  I4h  gave  much  hotter  results. 
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In  the  fitting  with  c  »  1.72,  there  were  still  two  somewhat  errant  orbits, 
at  MJD  42238  and  42324.  At  each  of  these  epochs  there  were  several  earlier  runs 
(nearly  as  good  as  the  final  orbit) ,  on  which  the  values  of  e  differed  by  up  to 
3  sd  from  the  value  on  the  final  orbit  used  to  obtain  '•  sd  of  e  was 

therefore  increased  by  a  factor  of  /To  on  these  two  orbits.  With  this  change, 
the  value  of  e  fell  from  1.72  to  1.08:  this  fitting  still  left  room  for 
improvement,  however,  because  the  largest  weighted  residual  (1.6)  was  for  the 
orbit  at  MJD  42245,  for  which  the  sd  had  already  been  increased  by  a  factor  of  10. 
In  least- squares  fitting,  the  largest  weighted  residual  has  the  greatest  effect 
on  the  sum  of  squares,  and  it  would  be  absurd  to  allow  this  ’discarded'  orbit  to 
exert  so  great  an  influence;  so  its  sd  was  increased  by  a  further  factor  of  1.6. 

The  resulting  (and  final)  fitting  of  the  values  of  e^  by  THROE,  shown  by 
the  unbroken  line  in  Fig  11,  gave  e  *  1.03  and  yielded  determinate  values  for 
all  four  of  the  lumped  harmonics,  as  follows: 


18.0  ±  10.0 

6- 1  » 1 

.0  s3,  - 

26.8  ±11.1  , 

(18) 

-  38.2  ±  4.2 

6^. 3 
io  s3!  - 

19.4  ±  2.8  . 

(19) 

In  Fig  11,  the  weighted  residuals  for  the  three  values  with  grossly  inflated 
standard  deviations  (shown  by  squares)  are  1.1,  1.0  and  0.6:  their  rms,  0.93,  is 
quite  close  to  the  overall  rms,  1.03.  So  these  three  orbits  do  not  significantly 
affect  either  the  fitting  or  the  sd  of  the  lumped  coefficients.  When  these 
three  orbits  were  omitted  and  the  remaining  14  orbits  were  fitted,  the  value  of 
e  increased  from  1.03  to  1.06,  and  the  values  of  the  lumped  coefficients  changed 
by  less  than  1 0%  of  their  sd  in  (18)  and  (19). 

The  effect  of  taking  the  point  M  of  maximum  density  at  16h  local  time 

instead  of  1 4h  was  tested  by  making  a  further  THROE  run  with  the  values  of  e2 
suitably  altered.  The  value  of  e  was  1.58,  as  compared  with  1.03  with  M  at 
14h.  This  comparison  is  not  entirely  fair  because  the  accuracies  of  some  orbits 
were  adjusted  to  improve  the  fitting  when  M  was  at  14h  local  time.  When  M  is 
at  1 6 h ,  the  orbit  at  MJL)  42324  no  longer  requires  its  sd  increasing  by  a  factor 

of  'To,  but  the  value  at  MJD  42238  requires  a  further  relaxation:  so  the  sd  of 

th  first  was  divided  by  /T5  and  that  of  the  second  multiplied  by  /l0.  A  new 
filting  by  THROE,  with  these  changes,  give  t  =  1.34  ,  still  considerably  higher. 
Thus  the  results  are  sensitive  to  the  assumed  local  time  of  maximum  density,  and 
a  better  fit  is  achieved  with  M  assumed  to  be  at  I4h. 
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6.5  Simultaneous  fitting  of  i  and  o 

The  SIMRES  computer  program'**  allows  a  simultaneous  fitting  of  i  and  e 
for  one  or  more  satellites.  If  the  e  terms  in  equation  (8)  are  important,  it 
is  clearly  preferable  to  make  a  simultaneous  fitting,  with 

h ,  q)  “  <1,  0),  (1,  i)  and  (1,  -1).  But  for  I973-27B  the  e  terms  in  equation 
(8)  are  not  expected  to  be  large,  so  the  advantage  of  a  simultaneous  fitting  is 
in  doubt  and  must  be  judged  empirically. 

The  fitting  wus  made  with  the  weighting  of  the  e  values  degraded  by  a 
factor  of  1.58,  which  is  the  ratio  of  the  values  of  c  on  the  fitting  of  e  und 
i  by  THROE  used  as  input  for  SIMRES.  The  values  of  the  lumped  coefficients  in 
the  simultaneous  fitting  were 


6-0.2 

6-0.2 

i°  c31 

-  1.46  *  0.19 

io  s3,  - 

-1.12 

±  0.25  , 

(20) 

6-1. » 

»o  c31 

-  20.8  ±  10.1 

6-1. 1 

io  s3!  - 

24.7  ± 

11.4 

(21) 

6-“  l .  3 

.o  c31 

-  -  36.1  ±  4.3 

6--1.3 
io  s31  - 

-H 

CO 

2.8 

(22) 

The  fittings  of  i  and  e  using  these  coefficients  are  shown  by  the  broken 
lines  in  Figs  8  and  11.  The  SIMRES  fitting  of  i  is  considerably  worse  than 
the  fitting  of  i  alone  with  (y,  q)  =  (1,  0),  the  sum  of  squares  of  weighted 
residuals  in  inclination  being  6.7  instead  of  4. 4.  With  e  the  simultaneous 
fitting  is  almost  as  good  as  the  fitting  of  e-  alone. 

6 , 6  Discussion 

It  is  difficult  to  decide  whether  to  recommend  the  values  from  the 

individual  fittings,  equations  (17)-(19),  or  the  values  from  the  simultaneous 

fitting,  equations  (20)— (22) ,  or  a  mixture  of  the  two.  If  the  results  from  i 

and  e  were  equally  reliable,  the  simultaneous  fitting  would  be  preferable.  But 

if  i  is  more  reliable  than  e  ,  the  results  obtained  from  f'tting  i  alone 

could  be  'contaminated'  by  being  forced  to  accept  less  reliable  values  of  the 

(> ,  q)  “  (I,  ±1)  terms  dictated  by  the  variations  in  e  .  Noting  that  (a)  the 

correction  to  e  was  very  sensitive  to  the  assumed  form  of  the  day-to- 

DN 

nivht  variation  in  air  density,  (b)  three  values  of  e  failed  to  fit,  and  (,c) 
tv.  >  others  were  given  increased  sd,  we  have  to  conclude  that  the  fitting  of  i 
is  more  reliable  than  that  of  e  .  The  probable  difference  in  reliability  seems 
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great  enough  to  justify  a  preference  for  the  fitting  of  i  alone,  equations  (17) 
This  opinion  -  and  it  is  only  an  opinion,  not  a  certainty  -  is  reinforced  by  the 
fact  that  the  unbroken  curve  in  Fig  8  fits  much  better  than  the  broken  curve, 
which  is  disturbingly  oscillatory.  For  eccentricity,  the  choice  is  not  important 
since  the  values  (18)  and  (19)  are  so  close  to  (21)  and  (22).  If  i  is  more 
reliable,  the  fitting  of  e  should  benefit  from  a  contribution  by  i  and  the 
values  from  the  simultaneous  fitting,  (21)  and  (22),  are  therefore  preferred. 

Of  the  three  pairs  of  lumped  coefficients,  (17),  (21)  and  (22),  the  third 
pair  may  seem  the  most  accurate  because  of  their  proportionally  smaller  sd.  But 
it  is  probably  fairer  to  compare  the  absolute  values  of  each  sd  after  division 
by  the  largest  numerical  coefficient  in  equations  ( 14)— (16)  respectively.  If 
equation  (14)  is  written  in  descending  numerical  order  of  coefficients,  we  have: 


C40, 3 1  “  0 ’ 799C42 , 3 1  "  °'701C38,31  + 


-0,2 
0.01 12C 


(9.0  ±  3.6)  x  10 


-0,2 

on  using  the  value  of  C3]  from  (17),  Similarly,  for  S  , 


>(23) 


S40,31  "  °*  799SA2 ,31 


-  0.701S 


38,3! 


(-  13.6  ±  2.2)  x  lo 


Equation  (15),  treated  in  the  same  way,  yields 


C41,3.  -  °-80lC43,31  -  °-707C39,31  +  —  “  <*  21 1  *  I0)  *  10 


S41  ,31  -  °-801S43,31  -  °-707S39,31  +  *"  =  C-  25  *  12)  -  10 


-9 


(24) 


Equation  (16),  with  similar  rearrangement,  gives 


C39  3|  -  0.893C41  31  -  0.617C37  +  ...  =  (134  ±  16)  *  10 


>  (25) 


S39 , 3 1  “  °-893sai  3 1  °-6  ’ 


...  =  (66  ±  10)  x  io 


-9 


Looked  at  in  this  light,  the  first  pair  of  coefficients  emerges  as  the  most 
accurate . 


The  Biopotential  coefficients  ol*  degree  i  are  expected  to  be  of  order  of 
-5  2  -9 

magnitude  !0  /(  ,  that  is,  of  order  6  *  10  for  (  '  4(1  .  So  the  lumped 

coefficients,  being  made  up  of  perhaps  10  or  15  significant  terms,  might  be 

-.0 

expected  to  be  of  magnitude  20  x  10  .  The  pairs  (23j  and  (24)  fulfil  this 

expectation;  but  the  pair  (25)  does  not,  and  the  C  coefficient  in  particular 
seems  unduly  large.  No  explanation  of  this  anomaly  can  be  offered  at  present. 

The  values  of  the  lumped  31st-order  coefficients  depend  critically  on 

inclination  and  the  numerical  factors  in  the  expressions  (I4)-(I6)  look  quite 

2 

different  to  those  in  the  corresponding  expressions'-  for  Proton  4  (inclination 
51.5°,  as  opposed  to  50.0°).  for  example,  C^j  is  multiplied  by  -  979.6  in 
(15)  and  by  -  534  for  Proton  4;  ^  is  multiplied  by  784.4  in  (15)  and  by 

256  for  Proton  4.  The  equation  from  Proton  4  corresponding  to  the  anomalous 
C  equation  in  (25)  is 

S).31  '  °-585V3,  -  °-7“S7,3l  *  •••  -  «■'  *  18 >  y  '0-9  ' 

So  Proton  4  also  gives  a  large  and  positive  value,  on  the  right-hand  side,  though 
the  numerical  factors  on  the  left-hand  side  differ  enough  to  rule  out  any  strict 
comparison. 

7  UPPER-ATMOSPHERE  ROTATION 

7 . 1  Perturbations  in  inclination 

The  values  of  inclination,  after  removal  of  relevant  perturbations,  should 
shoo  u  continual  decrease  produced  by  the  effects  of  upper-atmosphere  rotation, 
and  it  is  possible  to  determine  the  upper-atmosphere  rotation  rate  A  from  the 
variation  in  inclination,  Eig  6  shows  the  values  of  inclination  after  removal  of 
zonal  harmonic,  J,,  ,,  ,  and  lunisolar  perturbations,  as  described  in  section  6.3, 
with  the  curve,  of  Fig  8  drawn  in  across  the  resonance  region. 

In  assessing  the  significance  of  other  perturbations,  it  should  be 
remembered  (a)  that  only  one  of  the  62  orbits  has  a  nominal  accuracy  better  than 
30  in,  and  (b)  that  the  satellite  was  25  in  long,  and  that  the  observations  relied 
on  reflections,  either  radar  or  optical,  coming  from  any  part  of  it.  Because  of 
these  uncertainties,  perturbations  of  order  50  m  or  less  must  bo  regarded  as 
negligible.  Among  these  negligible  perturbations  arc  the  effects  of  Earth  and 
ocean  tides  (35  m) ,  solar  radiation  pressure  (10  m)  and  meridional  winds  (20  in) . 
The  unimportance  of  meridional  winds  may  seem  surprising,  because  such  winds  can 
have  appreciable  effects  on  an  eccentric,  orbit  when  the  perigee  is  at  particular 
local  times.  Meridional  winds  have  a  diurnal  variation,  however,  being  away  from 
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the  equator  by  day  and  towards  the  equator  by  night"  ,  and  for  a  nearly  circular 
orbit  like  that  of  1973-2711,  they  will  tend  to  cancel.  Even  if  the  meridional 
winds  are  consistently  in  one  direction,  however,  their  effects  still  tend  to 
cancel  on  a  near-c i rcular  orbit.  The  effect  of  a  constant  100  m/s  south-to-nor th 
wind  on  the  orbit  of  1973-27B  has  been  calculated:  the  maximum  change  in  inclina¬ 
tion  produced  was  0.00016°  (19  m) ,  which  must  be  regarded  as  an  upper  limit, 
since  a  constant  wind  is  unrealistic. 

Theoretical  curves  of  the  variation  of  inclination  with  time  for  1973-27B 

have  been  calculated  for  a  range  of  atmospheric  rotation  rates  A  using  the 

program  R0TATM,  which  is  based  on  a  lengthy  theoretical  equation  (quoted 

explicitly  by  Walker'^),  which  takes  account  of  atmospheric  oblateness,  and  the 

variation  of  scale  height  with  height.  The  only  other  effect  which  might  be 

...  27 

influential  is  that  of  the  day-to-night  variation  in  density,  but  theory 
indicates  that  this  effect  is  small  unless  the  orbit  is  close  to  certain 
conditions  of  resonance  between  the  motion  of  the  Sun  and  that  of  the  perigee: 
these  conditions  do  not  arise  for  I973-27B. 

7.2  Rotation  rate  from  June  1973  to  July  1974 

Between  launch  and  the  31:2  resonance,  the  theoretical  curve  with  A  =  1.09 
best  fits  the  points  in  Fig  6.  The  fit  is  entirely  satisfactory,  and  the  value 
of  A  is  quite  accurately  determined  by  the  end  points.  The  first  8  orbits 
serve  to  define  the  starting  value  as  50.0942°  ±  0.0002°.  The  final  value  can  be 
specified  in  two  ways:  either  (a)  from  the  group  of  three  accurate  independent 
values  near  MJD  42200,  which  gave  50.0390°  -  0.0002°  at  MJD  42210,  equivalent  to 
50.0385°  -  0.0002°  at  MJD  42232;  or  (b)  by  the  curve  of  Fig  8,  which  shows  that 
the  mean  of  the  oscillatory  fitted  curve  is  0.0006°  above  the  observational  value 
of  MJD  42232,  giving  50.0387°  ±  0.0002°  (on  doubling  the  nominal  sd) .  Both  (a) 
and  (b)  therefore  give  essentially  the  same  result  an!  if  we  assume  that 
i  **  50,0386°  ±  0.0002°  at  MJD  42232,  the  decrease  in  inclination  from  the  start 
is  0.0056°  ±  0.00028°,  giving  an  error  of  -5%  in  A  .  Thus  A  =  1.04  ±  0.05 
between  17  June  1973  and  4  July  1974.  This  value  is  averaged  over  local  time, 
because  (a)  the  orbit  is  so  nearly  circular  and  (b)  the  (ill-defined)  perigee 
moves  through  two  day-to-night  cycles  between  MJD  41S50  and  42232  (see  Fig  15). 

To  find  the  effective  height  y,  at  which  this  first  value  of  A  applies, 
we  note  that  the  mean  perigee  height  between  June  1973  and  July  1974  (see  Fig  lb) 

was  355  km,  where,  from  CInA  1972,  the  average  value  of  H  in  1973-4  was  50  km. 

20 

The  average  value  of  z  =  ae/H  is  therefore  0.7  and  theory"  indicates  that  v. 
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should  be  taken  at  Q.5H  above  perigee.  So  A  applies  at  a  mean  height  of  380  km, 
with  a  possible  error  of  about  5  km. 

The  curve  drawn  through  the  mean  values  of  A  obtained  in  previous 
2'  28 

analyses  ’  ,  reproduced  in  Fig  12,  shows  a  maximum  of  about  1.3  at  a  height  of 
350  km  with  a  rapid  decrease  to  about  1.0  at  a  height  of  400  km.  The  value 
obtained  from  197 3—27 B ,  namely  1,04  ±  0.05  at  380  ±  5  km,  is  below  the  curve,  and 
suggests  that  the  decrease  in  A  may  occur  at  a  height  about  20  km  lower  than 
the  curve  suggests  -  at  least  in  1973-1974. 

For  1 9 73—27 B ,  the  effective  height  y^  decreases  from  390  km  initially  to 
365  km  in  July  1974:  so,  if  the  curve  of  Fig  12  is  correct,  the  value  of  A 
should  increase  between  June  1973  and  July  1974.  A  slightly  lower  value  over 
the  first  six  months  and  a  slightly  higher  value  over  the  next  six  months,  would 
fit  the  values  in  Fig  6  just  as  well;  but  the  decrease  in  inclination  is  too 
small  to  justify  more  than  a  single  mean  value. 

7.3  The  31:2  resonance  (4  July  to  10  October  1974) 

In  the  analysis  of  the  resonant  variation  of  inclination  (section  6.3), 
better  results  were  obtained  with  A  =  1.3  than  with  A  =■  1.2  .  This  cannot  be 
called  an  'evaluation'  of  A  ,  but  does  show  that  the  observed  variations  arc 
consistent  with  a  value  of  A  of  1 ,3.  The  effective  height  y^  is  350  km,  and 
the  value  is  shown  by  a  diamond  in  Fig  12. 

7.4  Rotation  rate  from  10  October  1974  to  decay  (11  January  1975) 

The  value  of  inclination  at  the  end  of  resonance  is  of  course  just  as  well 
defined  as  that  at  the  beginning,  and  from  Fig  8  it  appears  that  the  last  value 
in  the  resonance  period,  at  MJD  42330  (10  October  1974),  coincides  with  the  mean 
of  the  oscillatory  fitted  curve.  So  this  value  of  inclination,  namely  50.0371° 
(after  removal  of  2  ’  luni-s°lar  anc*  zonal  harmonic  perturbations),  serves  as 
a  secure  starting  point  for  the  continuation  of  the  theoretical  curve.  The 
ruminal  accuracy  of  the  value,  0.0001°,  can  appropriately  be  doubled  to  allow  for 
other  sources  of  error,  such  as  the  averaging  of  Earth  and  ocean  tide 
perturbations . 

The  most  accurate  value  of  inclination  between  MJD  42330  and  decay  is  that 
at  MJD  42393,  namely  50.0294°  (after  removal  of  perturbations).  It  seems  appro¬ 
priate  to  select  A  so  that  the  theoretical  curve  passes  close  to  this  point:  if 
we  take  A  -  1.34  ,  all  the  observational  values  between  MJD  42330  and  MJD  42409 
are  satisfied  to  within  1.1  sd.  Fig  13  shows  the  observational  values  and  the 
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theoretical  curve  for  A  =  1.34  with  a  time  scale  wider  than  in  Fig  6.  Some 
of  the  values  beyond  MJi)  42409  are  also  satisfied  by  the  A  =  1.34  curve,  but 
the  last  four  values  are  far  above  this  curve,  and  it  is  obvious  that  a  lower 
value  of  A  is  required  al  the  end  of  Lho  life.  The  value  which  best  satisfies 
the  last  four  points  is  A  =  1.06  ,  and  this  curve  intersects  the  A  =  1.34 
curve  at  a  date  close  to  MJI)  42409  (1974  Dee  28).  Se  the  theoretical  curves 
chosen  to  fit  the  values  are  A  =  1.34  for  MJD  42330-42409  and  A  <*  1.06  for 

MJD  42409-42423.  The  change-over  or  'hinge-point*  at  MJD  42409  is  marked  with  a 

black  disc  in  Fig  13.  (Using  more  than  two  values  of  A  would  give  a  closer 

fit,  but  the  two-value  fitting  is  adequate  because  it  satisfies  all  the  points  to 

within  1.1  sd. ) 

The  accuracy  of  these  two  determinations  of  A  is  quite  difficult  to 
assess.  It  seems  reasonable  to  assume  that  the  curve  is  accurate  to  0.0003°  at 
M.TD  423'13,  where  the  observational  value  has  sd  0.0005°  and  several  other  points 
nearby  arc  adequately  Littod.  But  the  two  further  points,  aL  MJD  42402  and 
42409,  do  not  strengthen  the  accuracy:  they  merely  confirm  it.  So  the  A  =  1.34 
curve  will  be  taken  as  applying  between  MJD  42330  and  42393,  and  as  measuring  a 
decrease  of  0.0080°  in  inclination  wi'h  an  accuracy  of  0.00054°  (or  adding  in  the 
error  of  0.0002°  at  the  start).  This  gives  A  =  1.34  ±  0.09  .  Between  MJD  42330 
and  42393  the  perigee  height,  decreases  almost  linearly  from  316  km  to  206  km, 
with  u  mean  of  290  km  (see  Fig  16);  so  the  effective  height  y^  is  305  *  5  km, 
being  0.411  above  perigee,  where  H  r  40  km. 

The  A  -  1.06  curve  at  the  end  of  the  life  fits  the  last  four  points  very 
well,  and  since  their  mean  sd  is  0.0006°,  the  curve  should  be  accurate  to  about 
0.0003°  at  the  end  of  the  life.  In  assigning  an  error  to  the  starting  value,  at 
MJD  42409,  we  note  thal  if  the  A  =  1.34  curve  is  accurate  to  0.0005°  at 
MJD  42393,  tue  extra  error  Incurred  through  A  being  10%  in  error  thereafter 
would  be  0.0004°,  giving  a  root-square  total  of  0.0006°  at  MJD  42409.  If  we 
accept  this  assessment,  the  error  in  A  is  6%,  or  0.06.  It  could  be  argued 
that  the  task  of  assigning  a  realistic  error  is  impossible,  but  a  tentative 
estimate  is  still  helpful.  The  height  at  which  the  value  A  =  1.06  applies  is 
also  questionable:  the  mean  ot  the  first  and  last  perigee  heights  is  200  km,  but 
the  curve  is  largely  determined  by  the  last  four  values,  when  the  perigee  height, 
decreases  I  rum  201  to  157  km,  with  a  mean  of  180  km.  It  seems  reasonable  to 
take  the  mean  perigee  height  as  190  1  10  km.  An  addition  of  0.3H,  whore 

H  30  km  ,  gives  y  =  200  i  10  km. 
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7 . 5  Annlyais  in  terms  of  orbital  period 

When  e  <  0.006  and  z  <  0.8  ,  as  for  1973-278,  the  lengthy  theoretical 
equation*^  giving  the  change  Ai  in  inclination  when  the  orbital  period  Tj 
changes  by  AT^  may  be  simplified  to 

<0  * 

where  I^(z)  is  the  Bessel  function  of  the  first  kind  and  imaginary  argument,  of 

degree  n  and  argument  z  ;  the  factor  is  nearly  equal  to  I,  being  given  by 

*^F  ~  1  -  0.040A  for  197 3—  2 711 ;  and  c  expresses  the  effect  of  atmospheric 

oblateness,  with  c  ~  0.15  for  1973-27B,  Since  l.,/l^  •  0.07  for  z  *  0.8  , 

the  variation  of  i  with  T,  should  be  nearly  linear  if  A  is  constant;  thus, 

d 

when  the  variation  of  i  with  time  is  non-linear,  as  with  1973-2711  after 
resonance,  a  clearer  picture  emerges  if  i  is  plotted  against  ,  as  in 
Fig  14,  After  resonance,  z  0.4  ,  so  that  l.,/l.  <  0.02  ,  and  the 
oversimplification  of  assuming  that  constant  A  corresponds  to  constant  slope  is 
correct  to  *2".  The  unbroken  lines  in  Fig  14  correspond  to  the  A  -  1.34  and 
A  ”  1.06  curves  in  Fig  13,  and  both  have  been  continued  as  broi;«n  lines  beyond 
the  binge-point  at.  MJD  42409,  marked  by  the  black  disc.  The  need  for  two  slopes 
is  obvious,  and  although  the  exact  choice  of  slope  and  of  the  hinge-point  is 
arguable,  there  is  apparently  no  other  preferable  pair  of  hinged  lines  passing 
through  the  point  at  the  end  of  resonance. 

7 . 6  Discussion 

The  three  values  of  atmospheric  rotation  rate  trev/day)  obtained  from  this 
analysis  of  I973-27B  arc: 

A  »  1.04  t  0.05  at  height  380  A  5  Uni,  for  dune  lv,73  to  duly  l°74; 

A  -  1.34  t  0.09  at  height  305  i-  5  km,  for  10  October  to  28  Dec  l°74; 

A  -  1.06  t  0.06  at  height  200  1  10  km,  lor  28  Dec  l‘>74  to  II  dan  197S. 

All  three  values  are  averaged  over  local  time,  because  the  orbit  is  so  nearly 
circular.  tFor  the  third  value  the  local  lime  at  perigee  is  near  midnight,  as 
Fig  15  shows;  but  the  orbit  is  very  nearly  circular,  with  o  ■'  0.001  ,  and  the 
midnight  perigee  would  tend  to  equalize  the  drag  round  the  orbit..) 
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Tii o  values  of  A  obtained  hurt'  arc  not  especial ly  accurate,  because  the 
total  do ere a so  in  inclination  is  quite  small  (0.03°);  but  the  values  do  have  tho 
virtue  that  they  cover  a  wide  range  oi  height.  Tho  existing  curve  tor  tho 
variation  of  A  with  height,  reproduced  in  Fig  12  -  showing  A  near  1.0  at 
200  luu,  1.3  at  350  km  and  1.0  at  400  km  -  is  based  on  a  compilation  of  results 
from  many  different  satellites,  at  different  dates,  latitudes,  etc.  Now  this 
curve  has  lor  the  first  time  buen  tested  by  results  from  a  single  satellite,  with 
confirmation  of  its  general  form,  though  tho  decrease  in  A  at  heights  between 
350  and  400  km  may  have  occurred  at  a  lower  height  in  1973-4. 

Average  values  ot  A  lowur  than  those  given  by  the  curve  oi  Fig  12  have 
been  obtained  in  two  recent  papers'^’  11  for  heights  of  220-320  km  during  the 
years  1971-1973.  Those  results,  from  analysis  of  the  satellites  1 9 7 1  —  1 0 (■> A 
(inclination  6b°)  and  1971-1 813  (inclination  70°),  are  in  conflict  with  the  high 
value  ot  A  (1.34)  obtained  from  1973-2711  at  305  km  during  tho  IubI  three  months 
of  1974.  Tho  difference  may  be  a  latitude  effect,  resulting  from  the  lower 
inclination  of  1973-2711:  a  tendency  for  low-inclination  satellites  to  give  high 
values  of  A  has  been  noted  previously^1,  and  it  is  probable  that  most  of  the 
'super-rotation'  of  the  upper  atmosphere  occurs  at  latitudes  below  45°. 

7 . 7  Comments  on  the  accuracy  of  the  orbital  model  in  PROP 

There  is  one  further  question  which  arises  when  (as  hore)  the  valueH  of 
inclination  obtained  are  pressed  to  tho  limits  of  their  formal  accuracy  to 
extract  information  on  A  .  The  question  may  bo  posed  provocatively  in  the 
form:  "Is  the  sd  of  a  PROP  orbit,  at  n  particular  epoch  invalidated  by  the 
neglect  of  lunisolar  perturbations  within  the  l’KOl’  model  during  tho  5  days  or  so 
over  which  tho.  observations  extend?" 

Lunisolar  perturbations  for  1^7 3—2711  have  buen  computed  at  daily  intervals 

i  y 

using  PROD  .  For  eccentricity,  tho  results  show  that  lunisolar  perturbations 
are  negligible,  never  more  than  7  m:  this  is  because  o  itself  is  so  small.  For 
inclination,  the  variations  are  much  larger,  often  of  order  !>0 mover  live  days. 
The  orbits  likely  to  he  affected  are  the  most  accurate  ones,  namely  orbits  Nos, 

2,  8,  14,  15,  19,  23,  27  and  34.  for  which  Lite  sd  in  inclination  is  less  than 
50  m.  The  ntiuvfmwn  departures  of  i  from  its  value  at  epoch,  due  to  lunisolar 
perturbations  during  the  days  covered  by  tho  observations  for  these  nine  orbits 
are  2H  m,  89  m,  3  m,  74  m,  4.3  m,  45  m,  1(>  m,  19  m  and  38  m  respectively.  Since 
(a)  the  epoch  is  always  taken  as  near  as  possible  to  the  mid-point  in  lime  ol  the 
observation:-.,  and  (b)  the  perturbation  in  i  is  usually  mono  tonic  over  live 
days  (if  not,  it  is  snvll),  the  orbit  determined  will  tend  to  average  out  the 
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per turbat ions  not  included  in  the  I’KOl’  model,  so  that.  the  resultant  error  in 
inclination  at.  epoch,  dj,,  say,  should  not  be  more  than  about  j  ol  the  maximum 
perturbation,  as  recorded  above,  For  eight,  ot  the  nine  orbits,  the  value  of  d  , 

L.  u 

thus  defined,  is  less  than  ^  of  the  standard  deviation  o  of  the  inclination  as 
given  in  Table  1,  and  therefore  would  not  be  expected  to  affect  o  appreciably. 

The  one  doubtful  orbit  is  No. 8,  the  most  accurate  ol  all,  which 
unfortunately  also  has  the  largest  luuisolar  perturbation,  l'or  this  orbit  PROP 
gives  o  »  0.00012°  ,  which  is  equivalent  to  !4m,  and  the  luuisolar  perturbations 
in  i  at  00  hours  on  the  7  days  centred  on  epoch  (.expressed  as  equivalent 
distances)  are:  8'J  in,  bb  tu,  3.S  m,  0,  -3J  m,  -S1'  m,  —  (<(»  m.  The  variation  with 
time  is  nearly  linear,  except  lor  the  last  value,  so  the  averaging  effect  implicit 
in  the  orbit  determination  should  ensure  that  errors  1  roui  this  source  do  not  much 
exceed  10  m.  The  accuracy  can  also  bo  assessed  i  tom  a  different  argument:  it  is 
the  Ucwitt  camera  observations  that  are  responsible  for  o  being  as  small  as 
14  iu,  and  there  were  two  Nets  of  these  camera  observations,  at  2  hours  and  at 
27  hours  alter  epoch.  Since  those  observations  are  dominant,  the  orbit  will  be 
'most  correct'  about  14  hours  alter  epoch,  with  the  luuisolar  perturbations  well 
away  from  epoch  easily  falling  within  the  observational  accuracy  of  Llie  loss 
accqrato  observations  available  there.  (The  nominal  accuracy  of  the  observations 
is  5  m  l'or  the  Hewitt  camera,  and  200  m  or  more  lor  the  other  observations.) 

This  argument  indicates  the  importance  of  taking  the  epoch  lor  orbit  determination 
as  close  as  possible  to  the  Hewitt  camera  observations.  In  the  14h  between  epoch 
and  the  time  when  orbit  No. 8  is  'most  correct'  the  luuisolar  perturbation  is 
I1)  til:  to  be  realistic,  therefore,  the  standard  deviation  of  i  on  orbit  8 
probably  needs  to  be  doubled,  to  0.0002°.  No  change  in  the  analysis  is  required, 
since  0.0002°  has  already  been  taken  as  the  accuracy  tor  the  mean  of  the  tirst 
eight  orbits.  Nor  is  any  change  needed  in  Pig  <>,  where  even  a  doubled  sd  is 
still  Luo  small  to  show. 

8  ANALYSIS  OF  VARIATIONS  IN  KCCENTRICITY  HUH  TO  AIR  ORAL 
8 .  I  Decrease  in  a ( 1  -  o ) 

11  the  variations  in  eccentricity  due  to  zonal  harmonic  and  luuisolar 

perturbations  are  removed,  the  remaining  variation  in  e  (outside  the  resonance 

region)  should  be  caused  almost  entirely  by  the  client  of  atmospheric  drag.  It 

.  24  ... 

we  accept  the  theory  for  the  el  fort  ot  air  drag  ,  the  variation  m  e  can  be 
analysed  to  determine  either  (a)  values  ot  scale  height  11  ,  il  the  factor  F 
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defining  the  amplitude  of  the  day-to-night  density  variation*  is  assumed  known, 
or  (b)  values  of  F  ,  if  11  is  assumed  knuwn ,  Alternatively,  we  can  assume  that 
the  values  of  H  and  F  given  by  CIHA  are  correct  -  because  previous  orbit 

analyses *  indicate  that  the  CIHA  values  of  11  are  correct  to  ±10%  for 
1973-4  -  and  then  compare  the  variation  o£  eccentricity  predicted  by  theory  with 
the  observed  variation.  The  latter  approach  has  been  adopted  here,  because  the 
theory  for  an  atmosphere  with  day-to-night  variation  in  deusily  has  not 
previously  been  tested.  We  work  with  a ( 1  -  e)  rather  than  e  because  the 
zonal  harmonic  perturbation  can  be  more  accurately  removed  (the  third  harmonic 
perturbation  to  e  ,  Ae^  »  ±s  given  by  J^R  sin  i  sin  co/j^a  ,  so  that  aAe^ 
remains  of  constant  amplitude  as  the  orbit,  contracts). 

The  values  of  a(l  -  e)  from  the  62  orbits,  plotted  as  triangles  in 
Fig  15,  have  an  oscillation  caused  chiefly  by  the  effect  of  odd  zonal  harmonics 
in  the  geopotential.  The  values  of  a(l  -  e)  have  been  cleared  of  zonal  har¬ 
monic  and  lunisolar  perturbations  using  the  computer  program  PROD  with  integra¬ 
tion  at  1-day  intervals.  The  orbits  given  by  l’ROD  are  drag-free,  and  the 
integration  was  restarted  whenever  u  on  the  drag-free  orbit  departed  from  the 
real  value  by  more  than  3°.  As  explained  in  sect 'on  6.4,  the  perturbation 
removed  from  a(l  -  c)  was  the  value  indicated  by  PROD  at  the  correct  value  of 
to  rather  than  at  the  correct  epoch. 

The  values  of  a ( 1  -  e)  cleared  of  these  perturbations  are  denoted  by  q 
and  plotted  as  circles  in  Fig  15.  It  will  be  seen  that  Q  continually 
decreases:  as  expected,  its  slope  is  generally  least  when  the  local  time  at 
perigee  is  near  14h  (when  the  dens  icy  is  near  maximum)  and  greatest  when  the 
local  time  at  perigee  is  near  04h. 

24  .  . 

Theory  indicates  that  the  decrease  of  the  eccentricity  (cleared  of 

gravitational  perturbations)  should  be  governed  by  the  equation 


da 

dx 


(z)  ■*  I.(z)F  cos  <|> 

P  1 

^  Iq(z)  +  12  ^  f  ^  C0S  ^p 

1  +  0(e)} 

(27) 


where  if.  is  the  geocentric  angular  distance  between  the  perigee  and  the  point 
M  of  maximum  density  (as  in  section  6.4),  z  =  x/11  ,  and  x  is  the  value  of  ue 
cleared  of  perturbations,  so  that  x  =  a  -  Q  ,  whence  dQ/da  =  1  -  dx/da  .  Thus 
the  theory  gives 

*  The  definition  of  F  (see  equation  (A-l)  of  the  Appendix)  is  such  that  the 
ratio  of  maximum  to  minimum  density  is  (I  +  F)/(l  -  F) .  For  1973-2711  during 
most  of  its  life,  F  ~  0.5  ,  so  that  the  ratio  is  about  ). 
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da 


I,  +  5  (In  +  I->)F  cos  4> 

-1  1 . 0 .  2 _ 11  j| 

In  +  I,F  COS  Ip  1 

U  1  p 


+  0(0.005)}  , 


(28) 


on  dropping  the  argument  z  .  If  values  of  H  and  F  for  the  appropriate  values 
of  height,  solar  activity,  etc,  are  taken  from  C1RA  1972,  the  right-hand  side  of 
equation  (28)  gives  theoretical  values  of  dQ/da  ,  which  can  be  compared  with  the 
observed  values  of  AQ/Aa  on  orbits  taken  far  enough  apart  to  ensure  that  the 
change  AQ  in  Q  between  the  two  orbits  can  be  accurately  determined.  Although 
the  overage  sd  in  a(l  -  c)  is  80  ra,  the  correction  for  zonal  harmonic  and 
lunisolar  perturbations  may  introduce  errors  of  a  similar  order  because  of  the 
discrepancies  between  the  values  of  u  on  the  drag-free  and  actual  orbits. 

Thus,  when  two  values  of  Q  are  differenced  to  obtain  AQ  ,  an  error  of  about 
0.2  km  can  be  expected,  and  if  an  accuracy  better  than  5 3  in  AQ/Aa  is  required, 
values  of  AQ  of  A  km  or  more  should  be  used. 

The  histogram  in  Fig  17  shows  the  observational  values  of  AQ/Aa  , 
calculated  on  this  basis.  (The  first  value  of  AQ  is  only  2.6  km,  but  it  should 
still  be  accurate  to  ±53,  being  based  on  an  average  over  9  orbits.)  Fig  17  also 
shows  the  valuos  of  dQ/da  given  by  the  right-hand  side  of  equation  (28)  when 
values  of  H  and  F  are  taken  from  C1RA  1972  at  the  appropriate  satellite 
perigee  height,  as  given  by  Fig  16,  and  the  appropriate  exospheric  temperature, 
taken  constant  at  800K  (because  there  was  little  variation  in  solar  activity  in 
197 3—4 ) .  The  shape  of  the  curve  joining  the  calculated  values  of  dQ/du  is 
drawn  to  correspond  with  the  rather  erratic  variation  of  <|>  ,  shown  at  the  top 

of  Fig  17.  The  broken  curve  in  Fig  17  shows  the  values  of  dQ/da  given  by 
theory  with  no  day-to-night  variation  in  density  (F  -  0). 

Comparison  of  the  theoretical  and  observational  values  of  dQ/da  in 
Fig  17  shows  excellent  agreement,  usually  well  within  the  possible  ±53  error  in 
the  observational  values.  The  resonance  region  ought  to  be  excluded  from  the 
comparison,  because  other  perturbations  occur;  but  even  here  the  agreement  is 
quite  good,  except  on  the  first  value  of  AQ/Aa  (which  is  affected  by  the  large 
increase  in  e  between  MJD  A2232  and  42251  visible  in  Fig  11). 

The  theoretical  curve  with  F  -  0  ,  shown  by  the  broken  line  in  Fig  17, 
provides  a  good  average  value;  but  the  theory  with  day-to-night  variation  in 
density  is  obviously  touch  better.  Its  departure  from  the  mean  curve  is  approxi¬ 
mately  proportional  to  F  ,  so  (if  we  accept  the  theory)  we  can  ask  whether  F 
is  Loo  small  or  too  largo.  The  answer  seems  to  be  that  F  is  approximately 
correct,  thus  confirming  again  the  general  accuracy  of  CIRA  1972. 
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The  values  of  dQ/da  in  Fig  17  beloro  resonance  are  on  average  about  27. 
lower  than  the  observational  values  ot  AQ/Aa:  this  bias  could  be  corrected  by 
an  increase  of  about  4 Z  in  the  assumed  value  of  11  .  Although  the  observational 
values  have  a  57.  sd,  the  bias  ol  2%  does  persist  over  six  successive  values  ol 
AQ/Aa  ,  so  it  is  probably  significant.  if  so,  the  conclusion  is  either  that  the 
CIRA  values  of  11  are  about  4Z  too  low  or  that  11  should  be  evaluated  at  a 
height  yH  above  perigee  (where  y  ~  0.5  here).  The  latter  seems  more  likely  to 
be  the  correct  conclusion  and  it  fits  in  with  the  simplified  theory  of  the 
Appendix,  where  11  has  to  be  evaluated  at  the  mean  of  the  perigee  and  apogee 
heights.  This  would  imply  that  thu  treatment  of  11  in  an  atmosphere  where  H 
varies  with  height^  can  be  carried  over  into  the  theory  for  an  atmosphere  with 

day-to-night  density  variation.  This  latter  theory  has  so  far  only  been 
24 

developed  on  the  assumption  of  constant  11  . 


8. 2  Variation  of  perigee  height 

The  distance  of  perigee  from  the  Earth's  centre  is  7n( 1  -  e)  +  dr  >  ,  where 

29  ' 

dr  is  a  small  perturbation  which  may  he  written 
P 


dr 


1  ■  72foK 
a(l"-"u) 


2. 

sm  1 


-e(3-x  sin^i)  +  e4(4-5  sin4i)  -  2(1  -  e  +  e‘t)sin*'i  sin4w 


2.. 


+  0(e 


3>] 


km 


(29) 


For  1973-278,  with  error  less  than  30  m, 


dr  -  0.97  -  1.94  siu“w  km  .  (30) 

P 

The  perigee  height  y  is  obtained  by  subtracting  the  local  Earth  radius 

P  2  2 

R^(-  6378.14  -  21.38  sin  i  sin  ui)  from  thu  perigee  distance.  Thus,  tor  1973-2711, 

y  “  u( 1  -  e)  -  (6371.86  +  5.31  cos  2m)  .  (3!) 

P 

The  values  of  y  are  plotted  in  Fig  16,  which  bus  already  been  used  in  choosing 
values  of  11  and  in  assessing  the  effective  height  y^  at  which  A  applies. 

The  high  peaks  in  Fig  16  occur  when  perigee  is  at  maximum  latitude  south;  the 
lower  pi aks  at  northern  apex. 

8 . 3  Values  of  ore entricity  neat  decay 

The  values  of  eccentric  y  in  Fig  3,  alter  tracing  out  a  spiral  locus 
through  most  of  the  satellite's  life,  shilt  to  a  dilteront  path  during  i  ite  last 
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8  days,  because  (a)  the  decrease  in  e  due  to  drag  begins  to  overcome  the 
increase  due  to  odd  zonal  harmonics,  and  (b)  the  day-to-night  variation  in  air 
density  begins  to  have  an  increasing  influence  on  the  argument  of  purigee. 


The  general  equations  for  the  variation  of  the  eccentricity  of  a  near- 

circular  orbit  in  an  atmosphere  with  day-to-night  variation  in  density  involve 

lengthy  expressions  with  Fresnel  integrals  of  three  different  arguments.  The 

equations  simplify,  however,  as  decay  approaches  and  indicate  that  z  ■*  1 1'  cos  4> 
-1  o  1 

and  u>  -*■  tan  (-  cos  i  tan  M'),  where  M'  -  Si  -  L  -  30  ,  if  the  maximum  of  the 

daytime  bulge  is  at  I4h  local  time  and  on  the  equator.  Here  L  is  the  solar 

longitude,  which  increased  from  282°  on  3  January  1975  to  290°  on  11  January 

1975  (the  day  of  decay).  The  theoretical  equations  also  indicate  that 


cos  <f> 


2  2  2,  i 

(cos  i  +  sin  i  cos  «') 


at  decay. 


This  theory  has  not  been  subjected  to  a  practical  test,  because  no  high- 
drag  circular  orbit  has  been  determined  with  sufficient  accuracy  at  daily 
intervals.  It  is  therefore  useful  to  test  the  theory  by  trying  to  interpret  the 
variations  of  e  and  w  for  I 973-27B  in  its  last  days.  The  limiting  values 
quoted  in  the  previous  parug  aph  may  be  regarded  as  'end-values'  towards  which 
the  orbit  would  be  moving  in  the  absence  of  the  odd  zonal  harmonic  perturbations. 
Lot  the  suffix  F  denote  these  'final'  values,  so  that,  on  any  ol  the  last  days 
of  thu  satellite's  life,  wc  can  obtain  values  of  z^,,  u^,  and  from  the 

equations 


zK  -  |F  cos  <!>  1 

(32) 

■  tan  ' (-  cos  i  tan  M') 
r 

(33) 

,  2.  .  2.  2  , .  J 

i  ,,  -  (cos  i  +  siu  x  cos  M')' 

pi 

(34) 

Values  of  z,,  and  uu  have  been  calculated  for  each  of  the  last  7  orbits  and 
F  F 

are  plotted  as  triangles  linked  by  a  broken  line  in  Fig  18,  in  the 
(z  cos  w,  z  sin  m)  plane. 

The  corresponding  observational  values  of  z  cos  w  and  z  sin  w  from  the 
lust  7  orbits  (with  one  exception)  arc  indicated  by  circles  in  Fig  18.  The 
exception  is  7  January,  for  which  the  values  of  e  and  w  in  Table  1  give  a 
non-conforming  point,  indicated  by  a  cross;  this  point  could  be  omitted,  but 
instead  it  has  been  replaced  by  values  of  e  and  u  from  a  previous  orbit  at 
the  same  epoch,  for  which  10^e  -  1.248  and  m  -  31.5°  .  These  are  the  values 
indicated  by  the  circle  in  Fig  18.  (This  orbit  is  also  used  iu  Fig  3,  with  the 
orbit  of  Table  1  as  a  cross.) 
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In  the  absence  of  air  drag,  the  future  variation  of  z  ami  u>  in 
would  be  (approximately)  motion  round  a  circle  centred  at  the  point  [U, 

This  point  is  the  equivalent  of  the  point  C  in  l'ig  3,  but  it  is  no 

longer  fixed  because  11  decreases  steadily,  from  32  hat  on  3  January  to  19  km  on 

11  January.  The  single-headed  arrows  in  Fig  18  show  the  expected  progress  along 

the  appropriate  circle  in  one  day  (or  3  days  for  the  first  point)  on  the  basis  of 

a  6°  duily  advance.  Drag  has  the  effect  of  ' attracting'  the  point  away  from  this 
path  and  towards  the  relevant  triangle:  this  'attraction*  is  suggested  by  double¬ 
headed  arrows,  their  length  being  (arbitrarily)  increased  as  decay  approaches. 

At  first,  drag  bus  little  effect  and  the  actual  path  is  close  to  the 
circle;  but  drag  gradually  takes  command,  and  by  January  11.0  the  drag  and 
gravitational  effects  ure  in  opposite  directions,  with  the  drag  dominant.  As 
height  decreases,  F  (and  hence  z^,)  tends  to  zero  and  so  both  loci  in  Fig  18 
are  heading  towards  the  origin. 

Thus  the  actual  variations  of  z  and  u>  are  as  would  be  expected  on  the 
basis  of  the  theory. 

The  theory  also  predicts  that  a  nearly  circular  orbit  decaying  in  an 
utraosphorc  with  a  duy-Lo-night  variation  in  density  will  adjust  its  perigee 
position  so  that  the  perigee  moves  towards  the  point  of  minimum  density  in  the 

orbital  plane.  When  perigee  is  at  the  minimum-density  point,  always 

o  o  ^ 

exceeds  90  but  cannot  roach  180  unless  the  orbital  plane  happens  to  pass 

through  the  minimum-density  point. 


Fig  18 
0.8a  \ 
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In  Fig  19  the  triangles  show  the  values  ol  «|>  in  the  last  days  ol  the 

satellite's  life,  from  equation  (34),  which  assumes  that  the  point  of  maximum 

(or  minimum)  density  is  on  the  equator.  The  circles  in  Fig  19  show  values  of 

^  on  the  actual  orbits. 

P 


Fig  19  provides  strong  support  for  the  theory,  since  the  real  values  of 
ifi  are  not  only  i  90°,  as  expected,  hut  also  within  13°  ol  the  corresponding 
.  In  early  January  the  point  of  minimum  density  is  likely  not  to  be  on  the 
equator,  as  assumed,  hut  at  a  latitude  perhaps"*'  as  high  as  3l)°N;  the  effect  of 
such  a  change,  may  be  roughly  assessed  by  taking  a  smaller  value  of  inclination 
in  equation  (34).  If  i  is  taken  as  33°  instead  of  30l  ,  the  curve  of  d*  j.  in 
Fig  19  moves  upward  and  nearly  coincides  with  41  »  ‘ls  shown  by  the  dotted  lino. 
So  even  the  small  difference  between  $  and  d  ,,  in  Fig  19  can  he  explained. 

p  pi 
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9  CONCLUSIONS 
9  .  I  Orbit  determination 

Ail  average  accuracy  of  about  90  ui  radially  and  across-track  was  achieved  in 
62  orbit  determinations,  based  on  the  observations  listed  in  Table  2.  Accurate 
daily  orbits  were  determined  in  the  six  days  before  decay,  with  the  aid  of  NORAD 
and  US  Navy  observations. 

There  was  unacceptable  scatter  among  the  values  of  inclination  originally 
obtained.  This  scatter,  which  was  eventually  eliminated,  resulted  partly  from 
the  slight  overweighting  of  geographically  isolated  observations,  and  partly 
from  the  inherent  defects  of  least-squares  fitting. 

9 . 2  The  31:2  resonance 

The  variations  in  inclination  and  eccentricity  during  31:2  resonance 

(July-Octobor  1974)  were  successfully  analysed,  after  some  difficulties  over  the 

removal  of  perturbations  in  eccentricity  caused  by  the  day-to-night  variation  in 

density.  The  analysis  yielded  values  of  six  lumped  geopoLeutial  harmonics  of 

order  31.  In  these  lumped  coefficients,  harmonies  of  degree  up  to  60  are 

significant  and  the  lumped  values  should  provide  a  test  of  the  accuracy  of 

coefficients  of  order  31  and  degree  31  to  60  in  future  evaluations  of  the 

30 

geopotential  which  go  to  order  and  degree  60  or  more.  (Current  models  do  not 
go  beyond  order  and  degree  36.) 

9.3  Atmospheric  rotation 

Three  good  values  of  the  average  atmospheric  rotation  rate  have  been 
obtained  over  a  wide  range  of  height.  Those  show,  for  the  first  time  from  a 
single  satellite,  the  peak  in  the  rotation  rate  aL  heights  of  300-350  km  Going 
upwards  in  height,  the  values  are  1.06  i  0.06  rov/day  at  200  ±  10  km  in  January 
1975,  1.34  ±  0.09  rev/day  at  305  ±  5  km  for  October-December  1974,  and 
1.04  ±  0.05  at  380  ±  5  kin  for  June  1973-July  1974. 

9 . 4  Variations  in  eccentricity 

During  most  of  its  life  the  satellite  was  moving  through  an  atmosphere  with 

a  large  day-to-night  variation  in  density,  which  strongly  influences  the  orbital 

eccentricity :  the  day-time  maximum  density  exceeded  the  night-time  minimum  by  a 

factor  of  about  3,  and  the  ellipticity  of  the  orbit  was  nearly  equal  to  the 

ollipticity  of  this  day-to-night  'bulge'.  The  variation  in  a(l  -  e)  was 

2  4 

predicted  using  theory*'  previously  untested  and  values  of  atmospheric  parameters 
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from  C1RA  197”:  the  predicted  values  agreed  well  with  the  observed  variations  in 
a(l  -  e)  duo  to  drag,  thus  vindicating  both  the  theory  and  CUu\ 

The  daily  orbits  near  decay  were  uccurate  enough  to  provide  a  practical 

Q 

test  of  the  specialized  (and  also  untested)  theory  for  near-circular  orbits 
about  to  decay  in  an  atmosphere  with  day-to-night  variation  in  density.  The 
theory  gives  specific  end-values  for  o  and  w  and  predicts  that  perigee  should, 
in  defiance  of  the  normal  effects  of  the  harmonic,  move  towards  the  point  of 

minimum  density  on  thu  orbit.  Both  these  predictions  were  tested  and 
satisfactorily  confirmed. 
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CORRECTION  TO  ECCENTRICITY  TO  REMOVE  THE  EFFECT  OF  THE 
UAY-TO-NICHT  VARIATION  IN  AIR  DENSITY 

Quite  realistic  models  now  exist  which  specify  the  variation  of  upper- 
atmosphere  density  with  height,  local  time  and  latitude.  Here  we  require  the 
best  simple  approximation,  and  we  assume  (as  in  Ref  24)  that  the  air  density  p 
at  a  fixed  height  y  depends  only  on  the  angular  distance  <t>  from  the  point  M 
of  maximum  day-time  density ,  taken  to  be  on  the  equator  at  14h  local  time,  and  is 
given  by 

p(y)  ”  p(y)(l  +  F  cos  <(>}  ,  (A- 1) 

where  p  is  the  mean  density  (when  4>  ■  90°)  and  F  is  a  constant,  which  has  a 
value  near  0.5  in  the  atmosphere  sampled  by  1973-27B.  When  <(>  “  90°  (mean  day- 
to-night  conditions),  we  assume  that  the  variation  of  density  with  height  can  be 
approximated  by  an  exponential. 


p(y)  -  pp  exp I 


y  -  y„ 


(A-2) 


where  p  is  the  mean  density  at  the  perigee  height  y  ,  and  II  is  a  mean 
P  P 

density  scale  height  -  the  value  which  gives  the  best  fit  to  the  actual  density 

at  heights  between  perigee  and  apogee,  that  is,  for  y  £  y  4  y^  +  2ae  .  In 

practice,  of  course,  the  variation  of  p  with  $  is  not  exactly  symmetrical,  as 

in  equation  (A-l),  and  the  variation  of  n  with  height  is  not  exactly 

exponential. 


On  combining  equations  (A-l)  and  (A-2),  we  have 


p  (1  +  F  cos  )  exp 
P 


(A-3) 


for  heights  between  perigee  and  apogee.  Evaluating  (A-3)  first  at  apogee  and 
then  at  perigee,  and  dividing,  we  have 


p  1  +  F  coo  t 

a  a 


1  +  F  cos  i|> 


’(-  f )  ' 


where  p  ,  p  are  the  densities  at  apogee  and  perigee  respectively,  and 
a  P 

<jj  ,  the  corresponding  values  of  <(>  .  In  calculating  the  changes  in  e  due 

a  p  9  'j  _ 

to  air  drag,  the  theory  calls  for  a  value  of  11  appropriate  to  a  height 
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+  i'll  ,  where  y  is  a  function  of  z  »  ae/H  .  For  orbits  with  z  >  2, 

y  =  1.5  ;  hut  for  1973-27B  at  31:2  resonance,  z  is  near  0.4,  for  which  y  -  0.4 

(Fig  t>.  3  of  Ref  23).  Thus,  for  1 9 7 3—27 B ,  H  should  be  evaluated  at  a  height 

y  +  0.411  ,  which  happens  to  be  half-way  between  perigee  and  apogee  height,  since 
P  _ 

0.4H  ae  .  This  means  that  H  may  be  taken  equal  to  H  ,  which  has  already 

been  defined  as  the  appropriate  mean  value  of  density  scale  height  between 

perigee  and  apogee.  Thus  z  =■  ae/H  and  since  <f>  ”  p  +  180°  ,  equation  (A-4) 

a  P 

becomes 


P 

P 


P 


1  -  F  cos  p 

■;  -  -v. - -r^-  exp(-  2z) 

1  +  F  cos  p 

P 


(A-5) 


Over  a  time  interval  of  a  few  days,  when  p  can  be  regarded  as  fixed,  the 
simplest  way  of  approximating  to  an  atmosphere  with  a  day-to-night  density 
variation  is  to  define  an  'effective  scale  height'  and  a  corresponding 

z  -  ae/H.  ,  such  that  a  simple  exponential  model  with  scale  height  H£  will 
give  the  same  result  as  (A-5),  Thus  we  require 


or,  on  using  (A-5) , 


r  -  c*p(-  2ke} 

p 


z  +  ^  In 


1  +  F  cos  p 

_ E 

1  -  F  cos  * 

P 


(A-6) 


(A-7) 


So  the  appropriate  value  of  z^,  at  each  epoch  can  be  evaluated  if  F  and  $ 
are  known. 

If  the  orbit  is  symmetrical  relative  to  the  'day-tiiue  bulge'  of  density,  -ie 

if  p  =0  or  180°  ,  the  variation  of  density  round  the  orbit  should  be  well 
P 

approximated  by  a  spherically  symmetrical  atmosphere  with  scale  height  11^,  .  In 
asymmetrical  situations  the  approximation  is  not  so  good,  but  fortunately  the 
efiects  of  high  drag  on  one  side  of  perigee  and  low  drag  on  the  other  (at  the 
same  numerical  value  of  eccentric  anomaly  E)  combine  to  produce  the  same  effect 
as  their  mean,  because  e  depends  only  on  cos  E  and  not  on  sin  E  (see  Ref 
23,  equation  4.11).  So  the  use  of  llj,  ,  though  obviously  inexaut,  is  a 
basically  stable  approximation. 

For  1 9 7 3—  2 713  in  July-October  1974,  the  mean  perigee  height  was  near  330  km, 
and,  guessing  H  =  50  km  ,  we  have  to  evaluate  H  at  a  height  0.411  -  20  km 
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above  perigee,  thut  is  at  350  km.  The  mean  level  of  solar  10.7  cm  radiation  was 
-22  -2  -1 

40  '  10  W m  Hz  ,  so  with  an  allowance  ol  40  K  for  geomagnetic  disturbance, 
the  mean  duy/night  exospheric  temperature  for  use  in  the  COSPAR  International 
Reference  Atmosphere  1972  is  800  K,  from  CIRA  1972 ,  page  265.  The  CIRA  tables 
then  give  H  -  45  km  for  y  »  350  km  ,  and  also  show  that  the  ratio  of  maximum 
d ay- t i me  to  minimum  night-time  density  is  2.8,  giving  F  -  0.47.  Thus  a/H  ■  ISO 
and  z  «  ae/H  »  150e.  Hence  values  of  z,  can  be  obtained  from  equation  (A- 7) 

li* 

tor  each  of  the  17  orbits  in  the  resonance  region,  taking  F  “  0.47  and  values 
of  $  from  Fig  10, 

Within  THROE,  the  uir-drog  correction  to  e  between  orbit  n  at  epoch  t 
and  orbit  (n  ♦  1)  at  epoch  t  j  ,  is  taken  ns 

4M0 

AuT  "  3M“  ylU)  *  (Vl  ’  tn)  *  (A"a) 


where 


in.  defined  in  section  6,1, 


M,  is  the  value  at  epoch  t  , 

n 

.  31 

is  a  ratio  of  Bessel  functions,  approximated  by 


and 


y,  (z) 


y,u)  - 


240z  +  153z2  +  68z3  +  24z4 
480  +  306/  +  196z2  +  80z3  +  24z4 


(A-9) 


Equation  (A— 8 )  assumes  a  spherically  symmetrical  atmosphere:  to 
effect  of  the  day-to-night  variation  we  need  to  replace  z,  by 
Thus,  since  THROE  already  makes  the  correction  Ae.^,  ,  the  i ixtini 
required  to  take  account  of  the  day-to-night  variation  is  Ae^ 


allow  for  the 
z^  in  (A-8) . 
correction 
,  given  by 


Ae 


ON 


(A- 1 0) 


The  total  correction  at  the  nth  epoch,  o  ,  is  found  by  summing  the  values  of 

^°DN  between  epochs  1  and  2,  epochs  2  and  3,  ....  epochs  (n  -  I)  and  n  .  The 

values  of  Ae^.  for  use  in  (A-10)  were  derived  initially  from  a  previous  THROE 

run  and  improved  by  iteration,  only  one  iteration  being  needed.  The  values  of 

e^jg  are  shown  in  Fig  9b,  the  full  line  giving  the  values  when  M  is  taken  al 

1 4h  local  time,  as  indicated  by  I'lh'A  1972,  and  the  broken  line  when  M  is  taken 

2  5 

at  161i  local  Lime,  as  indicated  in  .lucrli  i  a  ’  s  1477  atmospheric  model  .  Changing 
M  to  1 0h  signil  icantly  changes  e  ^  and  runs idcrabl  worsens  the  tit  ol  the 
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theoretical  curve  to  the  values  of  e  (sue  section  6.4).  Thus  the  simple  model 

used  here  works  better  if  M  is  taken  at  14h.  According  to  the  recent  MS1S 
32 

atmospheric  model  ,  the  local  time  for  maximum  density  varies  during  the  course 
of  the  year  between  14h  and  16h. 
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31 st-order  harmonic  coefficients  in  the  guopotuntial;  thusu  will  provide  a  crucial 
test  of  the  accuracy  of  future  goopotoutial  mr.de  1  s . 

The  variation  in  inclination  before  and  alter  resonance  hus  beun  analysed  to 
determine  the  average  atmospheric  rotation  rate  A  (rov/day).  Results  are: 

A  -  1.C4  ±  0,05  at  height  380  kin  between  June  1973  and  June  1974;  A  ■  1.34  ±  0.09  ut 
height  305  km  for  Oetohor-Deccr.bor  1974;  and  A  -  1.06  ±  0,06  at  height  200  km  in 
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